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The  goal  of  this  thesis  was  to  develop  a  better  method 
for  Air  Force  use  in  evaluating  Fresnel  region  power  density 
hazards  associated  with  paraboloidal  antennas.  1  wanted  to 
not  only  improve  on  the  aperture  illumination  model,  but 
improve  on  the  entire  procedure  for  calculating  the  power 
density  figures.  The  resulting  computer  program,  POWERDEN , 
utilizes  the  one-parameter  circular  aperture  distribution 
developed  by  R.  C.  Hansen,  and  I  want  to  publicly 
acknowledge  his  work  here.  I’ve  attempted  to  make  the 
program  as  user-friendly  as  possible,  although  it  is 
designed  for  use  by  those  people  already  familiar  with 
paraboloidal  radiation  and  the  related  issues  of 
electromagnetic  compatibility.  Besides  absolute  power 
density  and  distance  figures,  the  output  in  both  the  axial 
and  off-axis  modes  is  also  given  in  relative  terms  so  that 
graphs  based  on  the  data  can  be  more  easily  constructed.  I 
hope  that  my  efforts  have  resulted  in  a  useful  product  which 
the  Air  Force  can  actually  use  to  help  protect  its  people 
and  equipment  from  hazardous  radiation. 
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Abstract 

r 

A  computer  program  is  presented  which  calculates  power 
density  in  the  Fresnel  region  of  circular  parabolic 
reflector  antennas.  The  aperture  illumination  model  is  the 
one-parameter  circular  distribution  developed  by  Hansen. 

The  program  is  applicable  to  the  analysis  of  electrically 
large,  center-fed  (or  Cassegrain)  paraboloids  with  linearly 
polarized  feeds. 

The  scalar  Kirchoff  diffraction  integral  is  solved 
numerically  by  Romberg  integration  for  points  both  on  and 
perpendicular  to  the  antenna  boresight.  Axial  results 
correspond  with  those  published  by  Hansen.  Off-axis  results 
cannot  be  directly  compared  to  any  others  obtained  with  this 
illumination  model,  but  they  are  consistent  with  what  is 
expected  in  the  Fresnel  region  where  a  quadratic  must  be 
added  to  the  linear  phase  term  of  the  integral  expression. 

Graphical  results  are  presented  for  uniform  illumination  and 
for  cases  where  the  first  sidelobe  ratio  is  20,  25,  30,  and 
35  dB.  ■.-:>] 


AN  IMPROVED  METHOD  FOR  CALCULATING  POWER  DENSITY  IN  THE 
FRESNEL  REGION  OF  CIRCULAR  PARABOLIC  REFLECTOR  ANTENNAS 


I .  Introduction 

Overview 

Circular  parabolic  reflector  antennas  (paraboloids)  are 
used  extensively  by  both  the  military  and  civilian  sectors 
for  microwave  communication,  radar  tracking,  and  deep-space 
applications.  Relative  design  simplicity  coupled  with  high 
power  gain  and  symmetric,  pencil-beam  radiation  patterns 
make  them  well  suited  for  these  tasks.  However,  the  same 
energy  focusing  characteristics  that  help  make  the 
paraboloid  so  useful  can  also  result  in  hazardous  power 
densities  in  their  vicinity  that  are  of  concern  because  of 
potential  biological  damage  such  as  cataracts,  interference 
with  pacemakers  or  other  electronic  equipment,  and  ignition 
of  nearby  explosive  material  (Hansen,  1976b;50;  Farrar  and 
Adams,  1980:134;  Department  of  the  Air  Force,  1981:  1-1; 
Michaelson,  1980:40-49). 

In  most  instances  the  hazardous  power  densities  occur 
within  the  antennas’  Fresnel  region  (radiating  near  field). 
In  this  region,  power  density  calculations  are  more 
difficult  than  those  associated  with  the  far  field  where 
radial  field  components  are  negligible,  and  all  rays  follow 
nearly  parallel  paths  to  the  observation  point.  Power 
density  in  the  far  field  varies  simply  with  the  inverse 


square  of  distance  according  to 

S  =  PG/4TR2  ( 1  ) 

where 

S  =  power  density 
G  =  far  field  power  gain 
R  =  far  field  range 

Use  of  Eq  (1)  in  the  Fresnel  region  gives  answers  that  are 
too  large,  since  gain  is  degraded  in  that  region  (Hansen, 
1976b: 52).  Rather  than  being  constant,  Fresnel  region  gain 
fluctuates  in  a  sinusoidal  manner  due  to  constructive  and 
destructive  interference  among  aperture  field  components 
with  different  amplitudes  and  phases  (Jull,  1981:46). 

The  method  used  now  by  the  Air  Force  for  determining 
these  hazardous  Fresnel  region  power  densities  is  outlined 
in  Chapter  6  of  Technical  Order  31Z-10-4  and  assumes  an 
aperture  illumination  that  is  rotationally  symmetric  in 
amplitude,  uniform  in  phase,  and  linearly  polarized.  The 
model  is  of  the  form 

G(x)  =  (l-xZ  )N  (2) 

where 

G(x)  =  normalized  aperture  distribution 

X  =  p/a,  normalized  aperture  radial  variable 

p  =  aperture  radial  variable 

a  =  aperture  radius 

N  =  exponent  controlling  rate  of  amplitude  taper 

and  is  deficient  in  the  following  respects: 

1.  It  results  in  zero  amplitude  at  the  aperture  edge. 
This  implies  either  an  infinite  aperture  or  the  perfect 


interception  of  all  energy  radiated  by  the  feed  horn. 

Neither  condition  occurs  in  practice  (Sciambi,  1965:79; 
Rabello  and  Nobili,  1973:677;  Hansen,  1976b:50);  instead, 
the  amplitude  tapers  off  to  some  minimum,  non-zero  value. 

2.  The  taper  control  exponent  is  restricted  to  the 
integer  values  of  0,  1,  2,  or  3.  There  is  no  physical 
reason  why  this  parameter  should  be  restricted  to  integer 
values  and,  for  highly  tapered  distributions,  it  can  exceed 
the  value  of  three  (Skolnik,  1970:9-20;  Silver,  1949:194; 
Farrar  and  Adams,  1980:135). 

3.  The  procedure  for  estimating  the  taper  control 
exponent  is  tied  to  the  beamwidth.  This  relationship  is 
ambiguous,  and  the  estimate  is  checked  by  calculating 
antenna  efficiency  using  a  gain  factor  corresponding  to  each 
of  the  four  allowable  integer  values  of  the  exponent. 
Succesive  values  of  the  gain  factor  are  plugged  into  the 
efficiency  equation  until  the  value  falls  anywhere  within 
the  range  0.5  to  0.9.  This  procedure  is  imprecise,  and  it 
is  possible  that  none  of  the  four  gain  factors  will  result 
in  an  acceptable  efficiency  figure.  Beamwidth  is  not  a  good 
indication  of  any  particular  illumination  because  it  is 
related  to  antenna  size  as  well  as  the  field  distribution 
across  the  aperture,  i.e.  even  if  the  illumination  remains 
the  same,  the  beamwidth  can  be  halved  simply  by  doubling  the 
antenna  diameter.  The  rate  of  amplitude  taper  is  actually 
related  to  the  level  of  the  first  sidelobe  of  the  far  field 


pattern  (Farrar  and  Adams,  1980:134;  Skolnik,  1970:9-25). 


Problem  Statement 

The  goal  of  this  thesis  is  the  development  of  a  better 
method  for  Air  Force  use  in  calculating  the  Fresnel  region 
power  densities  associated  with  paraboloidal  antennas.  The 
method  will  utilize  an  aperture  illumination  model  that 
corrects  the  three  deficiencies  discussed  above  and  will  be 
incorporated  into  a  computer  program.  It  will  directly 
calculate  axial  power  density  along  the  antenna  boresight  as 
well  as  power  density  off-axis  so  that  a  radiation  hazard 
radius  based  on  existing  standards  can  be  determined. 

Assumptions 

1.  The  paraboloid  is  electrically  large.  This  condition 
is  almost  always  met  when  using  practical  antennas  at 
microwave  frequencies  and  explains  why  paraboloidal  design 
is  based  on  high  frequency,  optical  methods  such  as 
geometric  optics  and  ray  tracing.  Methods  such  as  GTD  which 
account  for  diffraction  from  the  rim  are  required  only  for 
extreme  accuracy  or  for  wide  angles  encompassing  the  minor 
lobes  (Balanis,  1982:612-613).  The  method  developed  in  this 
report  is  based  on  a  minimum  aperture  diameter  of  lOX. 

2.  The  paraboloid  is  focused  at  infinity.  Far  field 
operation  is  usually  the  goal  so  that  this  requirement  is 
met  in  the  majority  of  cases.  Reflector  antennas  are 
sometimes  defocused  for  the  purpose  of  maximum  energy 


transfer  as  the  target  approaches  or  comes  within  the 
antenna’s  near  field  (Kay,  1960:587). 


3.  The  aperture  illumination  is  rotationally  symmetric 
in  amplitude.  This  is  an  ideal  condition,  but  designers 
usually  do  try  to  realize  symmetric  distributions  (Hansen, 
1964:  64;  Rudge  et  al . ,  1982:42)  so  that  the  subsequent 
analysis  will  be  simpler. 

4.  The  aperture  illumination  is  uniform  in  phase.  This 
is  also  an  ideal  condition  that  can  never  be  totally 
realized  because  of  unavoidable  imperfections  in  the 
paraboloid’s  construction.  Although  deviations  from 
constant  phase  might  be  planned  for  beam  steering  or  some 
other  specialized  purpose  (Jull,  1981:26),  they  are  usually 
unplanned  and  result  in  lower  gain,  higher  sidelobes,  and 
wider  beamwidth  (Skolnik,  1980:233). 

5.  The  electric  field  is  linearly  polarized.  The 
assumption  of  linear  polarization  along  with  the  high 
frequency  approximation  addressed  earlier  allows  calculation 
of  the  diffracted  field  along  the  antenna  boresight  without 
regard  to  the  field’s  vector  nature,  i.e.  it  becomes  a 
scalar  diffraction  problem  (Silver,  1949:164).  Although  the 
cross-polarized  field  components  do  not  cancel  completely  at 
points  located  off-axis,  their  effect  is  usually  neglected 
since  most  practical  paraboloids  focus  energy  in  a  very 
narrow  angular  region  around  the  axis  of  symmetry. 

6.  The  effects  of  aperture  blockage  are  neglected.  The 
center  feed  horn  (or  hyperbolic  subreflector  in  the 
Cassegrain  case)  and  its  related  supports  are  unavoidable 
obstacles  which  alter  the  radiation  pattern  of  paraboloidal 


antennas.  Although  the  effect  of  the  supports  has  proven 
difficult  to  analyse,  it  has  been  shown  by  superposition  of 
aperture  illuminations  that  the  effect  of  the  central 
blockage  is  to  decrease  gain,  increase  the  level  of  the  odd 
numbered  sidelobes,  and  decrease  the  level  of  the  even 
numbered  sidelobes  (Clarke  and  Brown,  1980:207-209;  Rudge  et 
al.,  1982:142-145). 

7.  The  antenna’s  Fresnel  region  will  be  assumed  to  begin 
a  few  wavelengths  in  front  of  the  antenna  and  continue  out 
to  the  conventionally  agreed  upon  beginning  of  the  far  field 
at  a  distance  of  2D^/A,  where  D  is  the  paraboloid’s  diameter 
and  X  is  the  wavelength.  Out  to  a  few  wavelengths  away  the 
reactive  field  dominates  and  no  approximations  can  be  made 
to  the  field  integral  (Silver,  1949:170).  The  20* /X  limit 
is  based  on  the  criteria  that  no  phase  difference  on  axis 
can  exceed  jr/8  radians  or  22.5  degrees  (Skolnik,  1970:9-9). 
Although  it  may  be  inadequate  for  highly  tapered 
distributions  which  result  in  extremely  low  sidelobe  levels 
(Hansen,  1964:31),  this  far  field  criterion  is  widely 
accepted . 

8.  The  reader  is  expected  to  have  an  understanding  of 
basic  antenna  concepts  such  as  gain,  beamwidth,  aperture 
illumination,  sidelobe  level,  and  far  field  radiation 


characteristics. 


Scope 

The  method  developed  for  this  thesis  will  be  applicable 
to  standard,  center-fed  paraboloids  or  to  dual  reflectors  of 
the  Cassegrain  type  that  employ  conventional  pyramidal  or 
conical  feed  horns.  It  will  not,  in  general,  be  applicable 
to  systems  with  less  common  features  such  as  offset  or 
corrugated  feed  horns,  since  the  previous  assumptions  of  a 
rotationally  symmetric  aperture  distribution  and  of 
negligible  cross-polarized  components  along  the  boresight 
might  not  be  valid.  The  resulting  computer  program  is 
expected  to  be  a  useful  tool  for  Air  Force  use  in 
determining  the  location  of  hazardous  Fresnel  region  power 
densities . 

Method 

The  approach  taken  will  be  the  development  of  a  field 
integral  equation  that  is  valid  in  the  Fresnel  Region.  The 
aperture  illumination  model  will  be  the  one-parameter 
circular  distribution  developed  by  Hansen  and  discussed 
further  in  Chapter  3.  The  coupling  of  a  non-zero  amplitude 
pedestal  at  the  aperture  edge  and  a  symmetric  amplitude 
taper  makes  it  a  more  realistic  model  (Sciambi,  1965:79; 
Rabello  and  Nobili,  1973:677;  Hansen,  1976b:50,52;  Hansen, 
1964:50,64)  than  the  present  one  given  by  Eq  (2).  The 
computer  program  will  perform  the  necessary  numerical 
integration  using  the  Romberg  adaptive  trapezoidal  technique 
and  will  directly  calculate  power  density  at  discreet  points 
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both  on  and  off-axis.  Results  will  be  compared  and 
contrasted  with  the  theoretical  work  of  Hansen  and  others. 


Development 

Information  on  concepts  relevant  to  paraboloids  and  the 
Fresnel  region  in  general  is  given  in  Chapter  2  while 
Chapter  3  presents  the  results  of  a  literature  search  on  the 
subject  of  determining  power  density  in  the  Fresnel  region. 
Chapter  4  explains  the  design  and  use  of  POWERDEN ,  the 
computer  program  produced  as  a  result  of  this  thesis. 
Graphical  results  of  the  program’s  use  are  presented  in 
Chapter  5  where  they  are  analysed  in  relation  to  the  findings 
of  Hansen  and  others.  Chapter  6  presents  some  final 
thoughts  and  recommendations  for  further  study. 


Overview 

This  chapter  attempts  to  summarize  some  basic  concepts 
upon  which  the  determination  of  Fresnel  region  power  density 
is  based.  Paraboloidal  design  is  discussed  and  the 
characteristics  of  the  field  within  the  Fresnel  region  are 
expanded  upon  since  they  may  not  be  as  familiar  as  those 
which  apply  to  the  far  field.  The  concept  of  Fresnel  zones 
is  introduced  because  it  is  a  convenient  way  to  visualize 
the  process  which  results  in  the  sinusoidal  variation  of 
axial  power  density  in  the  Fresnel  region.  The  chapter 
concludes  with  a  discussion  of  the  approximations  which  are 
usually  made  to  the  field  integral  equation  when  the 
observation  point  is  within  the  Fresnel  region. 

Paraboloidal  Design 

The  principle  upon  which  the  design  of  paraboloids  of 

typical  dimensions  is  based  has  been  summarized  as  follows: 

It  has  been  shown  by  geometrical  optics  that  if  a 
beam  of  parallel  rays  is  incident  upon  a  reflector 
whose  geometrical  shape  is  a  parabola,  the  radia¬ 
tion  will  converge  (focus)  at  a  spot  which  is  known 
as  the  focal  point.  In  the  same  manner,  if  a  point 
source  is  placed  at  the  focal  point,  the  rays  re¬ 
flected  by  a  parabolic  reflector  will  emerge  as  a 
parallel  beam.  (Balanis,  1982:604) 

This  principle  is  illustrated  in  Figure  1  for  the  case  of  a 

circular  parabolic  reflector.  In  the  limit  of  zero 

wavelength  and  ideal  conditions,  a  spherical  wave  emanating 


from  the  focal  point  is  tranforaed  into  a  plane  wave  after 
reflection.  The  reverse  is  true  for  the  receive  mode  of 


Figure  1.  Operating  Principle  of  a  Paraboloidal  Antenna 
(Skolnik,  1970:10-3) 

operation.  These  excellent  energy  focusing  and  collimating 
prop>erties  are  part  of  the  reason  why  the  paraboloid  enjoys 
such  widespread  use. 

Another  reason  for  their  popularity  is  relative  design 
simplicity  which,  as  stated  above,  is  based  upon  geometric 
optics  or  ray  tracing.  Geometric  optics  is  a  high  frequency 
approximation  to  Huygen’s  principle  of  physical  optics; 
Huygen’s  principle  states  that  "each  point  on  a  primary 
wavefront  can  be  considered  to  be  a  new  source  of  a 
secondary  spherical  wave  and  that  a  secondary  wavefront  can 
be  constructed  as  the  envelope  of  these  secondary  spherical 
waves".  In  contrast,  geometric  optics  considers  each  point 
on  the  secondary  wavefront  to  have  a  direct  correspondence 
with  one  other  point  on  the  primary  wavefront  (Kraus, 
1984:525,529).  These  ideas  are  illustrated  in  Figure  2. 
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Figure  2.  Secondary  Wavefront  Formation  for  Physical  Optics 
and  Geometric  Optics  (Kraus,  1984:525,530) 

When  applied  to  the  analysis  of  paraboloids,  which  are 
typically  many  wavelengths  across,  geometrical  optics  can  be 
used  in  conjunction  with  the  aperture  integration  method  to 
find  the  radiated  field.  Ray  tracing  is  used  to  determine 
the  distribution  of  the  tangential  field  component  over  a 
plane  that  is  normal  to  the  axis  of  symmetry  and  that 
usually  passes  through  the  focal  point.  This  surface  is 
referred  to  as  the  aperture  plane,  and  the  field  is  assumed 
to  be  zero  outside  the  projected  area  of  the  reflector  onto 
this  surface  (Balanis,  1982:612).  Integration  over  this 
surface  and  aperture  distribution  is  then  performed  to  find 
the  radiated  field.  All  ray  paths  from  the  focal  point  to 
the  aperture  plane  are  of  equal  length  so,  under  ideal 


1  1 


conditions,  the  aperture  illunination  will  be  uniform  in 
phase.  This  is  not  the  case,  however,  with  the  amplitude 
portion  of  the  distribution.  The  power  density  of  spherical 
waves  from  the  focal  point  falls  off  with  the  square  of 
distance  while  plane  waves  present  after  reflection  suffer 
no  such  spreading  loss.  Thus,  components  reflected  near  the 
center  undergo  leas  attenuation  than  those  reflected  nearer 
the  edge,  and  a  natural  amplitude  taper  appears  in  the 
aperture  distribution  (Jull,  1981:66;  Clarke  and  Brown, 
1980:191).  Careful  design  of  the  feed  and  shaping  of  the 
reflector  surface  can  be  used  to  increase  or  decrease  this 
natural  taper  depending  on  the  desired  radiation  pattern. 

Fresnel  Region  Field  Characteristics 

The  Fresnel  region  extends  from  a  few  wavelengths  in 
front  of  the  aperture  plane  out  to  the  beginning  of  the  far 
field  at  a  distance  of  2D2/X.  In  the  far  field  any  point  is 
essentially  the  same  distance  from  all  the  points  on  the 
aperture  plane;  the  radiation  pattern  is  independent  of 
range  and  power  density  varies  simply  with  the  inverse 
square  of  distance  ( Eq  (1)).  The  Fresnel  region  "is 
characterized  by  the  onset  of  diffusion  outside  the 
boundaries  defined  by  the  extension  of  the  rays  through  the 
aperture"  (Silver,  1949:172).  All  ray  paths  to  a  field 
point  are  no  longer  parallel  so  that  constructive  and 
destructive  interference  effects  become  important.  The 


radiation  pattern  becomes  dependent  on  range  as  well  as 


angular  displacenent  with  typical  results  being  broadening 
of  the  main  beam,  blurring  of  sharp  nulls,  and  an  increase 
in  sidelobe  levels.  These  effects  of  finite  range  become 
more  pronounced  as  range  decreases.  They  can  be  accounted 
for  mathematically  by  the  addition  of  a  quadratic  phase  term 
to  the  linear  one  present  in  the  far  field  equation  (Jull, 
1981:31,34).  Axial  power  density  no  longer  varies  with  the 
inverse  square  of  distance  in  this  region,  but  fluctuates 
sinusoidally  at  an  increasing  rate  as  the  range  is 
diminished.  Gain  is,  in  general,  degraded  with  aperture 
distributions  that  display  the  highest  far  field  gain 
undergoing  the  most  significant  reduction  (Hansen,  1976b:5l; 
Jull ,  1981:51). 

Fresnel  Zones 

The  division  of  a  paraboloid’s  circular  aperture  plane 
into  Fresnel  zones  (Figure  3)  aids  in  understanding  why  the 
axial  power  density  fluctuates  sinusoidally  in  this  region. 
When  the  observation  point  is  in  the  far  field  it  is  nearly 


Figure  3.  Fresnel  Zones  (Clarke  and  Brown,  1980:154) 


the  same  distance  from  all  points  on  the  aperture  plane, 
and  there  is  no  interference  among  aperture  field 
components  that  are  out  of  phase.  As  range  decreases  a 
point  is  reached  where  the  path  length  difference  between 
the  aperture’s  center  and  edge  is  exactly  x/2.  At  this 
point  the  first  Fresnel  zone  encompasses  the  entire 
aperture,  and  axial  power  density  reaches  its  first  maximum. 
As  range  decreases  further,  more  and  more  zones  are  formed. 
The  path  length  difference  between  a  point  on  the  edge  of 
any  zone  and  the  aperture  center  is  an  integral  number  of 
hal f-wavelengths .  Thus,  for  every  contribution  from  one 
surface  element  in  a  particular  zone,  there  is  another  one 
from  an  identical  surface  element  in  the  next  zone  that  is 
180°  out  of  phase.  The  result,  as  illustrated  in  Figure  4 
for  the  case  of  uniform  illumination,  is  alternating  maxima 
and  minima  corresponding  to,  respectively,  the  formation  of 
an  odd  or  even  number  of  Fresnel  zones  (Silver,  1949:196- 
197),  (Note:  This  report  will  adopt  the  usual  convention  of 
presenting  all  plots  of  axial  power  density  vs  range  with 
values  of  both  normalized  with  respect  to  those  at  the  far 
field  distance  of  2D^/A  so  that  A  =  R/{2D^/A  )  in  Figure  4.  ) 
There  will  be  complete  cancellation  at  certain  axial  ranges 
only  for  the  particular  case  where  the  aperture  distribution 


is  uniform. 


If  the  antenna’s  aperture  distribution 
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Figure  4.  Power  Density  vs  Range  for  Uniform  Illumination 
(Jull,  1981:47) 


1949:170).  It  asumes  high  frequency  and  linear 


polarization . 


If  one  also  assumes  uniform  phase  across  the 


^  ^  V  -■  V  >  J* 


aperture  distribution,  then  this  equation  is  given  by 


na  exp(-jkr) 

G(p)  - 

4Trr 


[ ( jk+l/r )cos( n, r ) + jk]  p  dp  d0  ( 


where 


Ep 

G(p) 

k 

(n,r] 


scalar  field  at  observation  point 
aperture  distribution 
propagation  constant,  2ir/X 

distance  from  aperture  point  to  field  point 
angle  between  normal  to  aperture  plane  and 
the  r  direction 


Figure  5.  Aperture  Geometry 


The  traditional  approximations  for  the  Fresnel  region  are 
the  following  (Skolnik,  1970:9-8): 

1.  The  1/r  term  inside  the  second  set  of  brackets  is 
negligible  in  comparison  to  k. 

2.  The  1/r  term  in  the  first  set  of  brackets  is  set 
equal  to  the  constant  1/R  where  R  is  the  distance  from  the 
origin  to  the  observation  point. 

3.  The  angle  (n,r)  is  replaced  by  (n,R)  =  8  . 


4.  The  r  in  the  phase  term  exp(-jkr)  is  expanded  in  a 
power  series  and  all  terms  higher  than  the  second  order  are 
truncated  so  that  (Hansen,  1964:26) 

r  =  R  -  psin(  6  )cos  (4> -0  ) 

+  [p2/2R][l  -  sin*  (  #  )cos2  (0 -4> )  ]  (4) 


5.  If  field  points  at  wide  angles  are  not  considered 
then  the  Fresnel  Small  Angle  ( FSA )  approximation  neglects 
the  sin*(fl)  term  in  Eq  (4)  so  that  r  in  the  phase  term  is 
usually  approximated  by  (Hu,  1960:345) 

r  =  R  -  /)sin(  )cos(  0  )  +  p* /2R  (5) 


With  these  approximations  the  Kirchoff  integral  can  be 


written  as 


/  / 

exp[  jkpsin(  S  )cos(<P-^  )  Ipdpdrn  (6) 


whe  re 


Ep  =  8(6) 
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8(6)  =  (l  +  cos(6)]  [jk/4irR]  exp(-jkR) 


Eq  (6)  is  the  basis  for  the  method  developed  in  this  report. 
This  equation  differs  from  the  standard  far  field  equation 
due  to:  (1)  the  addition  of  a  quadratic  phase  term  so  that 
the  7r/8  phase  criterion  can  be  maintained  and  (2)  the  fact 
that  the  [l+cos(6)]  term  is  not  set  equal  to  two. 

Range  o f  Accuracy .  The  Fresnel  approximations  will 


give  accurate  results  for  observation  points  which  are  not 
too  close  to  the  aperture  plane  and  which  are  not  too  far 
off  the  antenna  boresight.  Close  to  the  antenna  the 


approximations  might  fail  due  to  consideration  of  phase 
errors  or  proximity  to  the  reactive  field  at  extremely  close 
ranges  (Hansen,  1964:30).  However,  only  for  extremely  small 
antennas  with  diameters  less  than  about  eight  wavelengths 
across  does  the  reactive  field  consideration  dominate 
(Skolnik,  1970:9-11).  For  axial  displacements  the 
approximations  have  been  shown  to  give  extremely  accurate 
results  as  close  as  two  diameters  away  from  a  uniformly 
illuminated  paraboloid  with  a  100  wavelength  diameter 
{ Rahmat-Sami i  et  al . ,  1981:162).  More  generally,  the  axial 
range  restriction  for  a  maximum  phase  error  of  T/8  has  been 
given  as 

R  >  0.5D(D/\)W3  (7  ) 


which  results  in  a  minimum  range  of  1.1  diameters  for  a  ten 
wavelength  aperture  and  2.3  diameters  for  a  100  wavelength 
aperture  (Hansen,  1964:30;  Skolnik,  1970:9-10).  Eq  (7)  is 
the  criterion  for  minimum  range  which  is  used  in  the  method 
developed  in  this  report. 

At  wide  angles  the  cross-polarized  field  components 
cannot  be  neglected,  and  the  FSA  approximation  mentioned 
earlier  ( Eq  (5))  is  no  longer  valid.  However,  since  most  of 
the  energy  is  concentrated  in  a  small  angular  region  by 
practical  paraboloids,  the  applicable  angular  range  is 
usually  sufficient.  Even  for  antennas  as  small  as  ten 
wavelengths  across,  the  FSA  approximation  gives  accurate 
results  for  angles  "covering  up  to  four  sidelobes  in  the  far 


field.  Since  this  approximation  is  absolute  with  regard  to 


angular  range,  the  approximation  improves  greatly  for  large 
D/X  with  regard  to  the  regions  of  greatest  interest  near  the 
principal  lobe"  { Rahmat-Sami i  et  al . ,  1981:162).  The 
general  relationship  between  antenna  dimensions,  axial 
range,  and  maximum  allowable  angle  off-axis  is  given  by 
(Zucker,  1966:685) 


Tra^ 

-  sin2(®)<<l  (8) 

Xz 

where 

a  =  antenna  radius 
X  =  wavelength 
z  =  axial  range 
$  =  angle  off-axis 

Eq  (8)  is  always  satisfied  in  the  far  field,  but  it  requires 
that,  in  the  Fresnel  region,  the  observation  point  be 
restricted  to  smaller  angles  as  the  range  decreases.  If  the 
relationship  is  made  more  specific  so  that 


-  sin* ( e)  <  0.08  (9) 

Xz 


then  a  maximum  off-axis  angle  of  26.8°  is  allowed  at  the 
beginning  of  the  far  field  where  z  =  2D*/X.  In  accordance 
with  the  findings  of  Rahmat-Samii  et  al .  cited  above,  this 
allows  coverage  of  almost  four  sidelobes  of  a  ten  wavelength 
antenna  at  that  distance.  Eq  (9)  is  the  criterion  for 
maximum  angular  displacement  used  in  the  method  developed  in 
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this  report. 
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III.  Literature  Review 


Overview 

This  chapter  presents  the  results  of  a  literature  search 
on  the  subject  of  determining  Fresnel  region  power  density. 
The  usual  procedure  followed  by  most  investigators  is  to 
solve  the  scalar  Kirchoff  diffraction  integral  after  the 
Fresnel  approximations  have  been  made  (Eq  (6)).  Constant 
phase  aperture  illuminations  that  are  rotationally 
symmetric  in  amplitude  are  usually  assumed  so  that  analysis 
is  simplified  and  numerical  integration  can  be  performed 
with  respect  to  only  one  variable.  The  method  of 
presentation  here  is  to  show  results  that  have  been  obtained 
assuming  progressively  more  complex  but  realistic  aperture 
distributions . 


Uni  form  1 1 lumination 

Aperture  distributions  are  usually  chosen  for  the  ease 
with  which  the  resulting  integration  can  be  performed  as 
well  as  for  their  similarity  to  practical  distributions.  A 
uniform  amplitude  distribution  is  not  very  realistic,  but 
has  been  investigated  extensively  because  it  results  in  a 
field  expression  which  can  be  evaluated  in  closed  form. 

This  illumination  utilizes  the  paraboloid’s  physical 
aperture  with  the  highest  efficiency,  and  other 
distributions  are  usually  evaluated  relative  to  this  case. 
In  the  far  field  it  results  in  the  highest  gain  and 


narrowest  beamwidth  of  any  constant  phase  distribution  but 
also  results  in  the  highest  sidelobes,  17,57  dB. 

The  Fresnel  region  gain  reduction  on  axis  for  a  uniform 
amplitude  distribution  has  been  solved  for  analytically 
(Jull,  1981:47;  Silver,  1949:199)  and  is  given  by 

G/Go  =  [3in(t)/t]2  (10) 

where 

G  =  Fresnel  region  gain 
Go  =  far  field  gain  at  infinity 
t  =  ka^ /4  z 

k  =  propagation  constant 
a  =  aperture  radius 
z  =  axial  range 

The  graphical  results  of  this  expression  have  already  been 
presented  in  Figure  4.  Using  Eq  (10),  the  power  gain  at  the 
beginnning  of  the  far  field  (z  =  20^ /X )  is  0.99Go  where 
Go  is  the  power  gain  at  infinity. 

Ampl i tude  Tape r  Without  Pedestal 

This  is  the  type  of  distribution  ( Eq  (2))  utilized  in  the 
method  currently  used  by  the  Air  Force  for  calculating 
Fresnel  region  power  density.  The  amplitude  taper  makes  it 
more  representative  of  practical  distributions  than  the 
uniform  illumination  just  discussed.  However,  as  mentioned 
in  Chapter  1,  zero  amplitude  at  the  aperture  edge  does  not 
usually  happen.  Nonetheless,  this  particular  distribution 
has  been  widely  used  as  a  model  for  paraboloidal 


i 1 lumination . 


A  computer  program  has  been  written  (Farrar  and  Adams, 
1980:134-137)  which  uses  this  illumination  to  calculate 
Fresnel  region  power  density.  In  contrast  to  the  present 
Air  Force  method,  this  program  does  not  restrict  values  of 
the  taper  control  exponent  to  integer  values,  and  they  are 
estimated  from  the  sidelobe  ratio.  Figure  6  shows  the 
fluctuations  in  Fresnel  region  power  density  assuming  this 
illumination  and  a  25  dB  sidelobe  ratio  (Hansen,  1964:38). 
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Figure  6.  Power  Density  vs  Range  for  Tapered  Illumination 
Without  Pedestal  (Hansen,  1964:38) 


Note  the  deviation  from  the  simple  monotonic  increase  as  the 
Fresnel  region  is  entered  at  ^  =  1.  Comparison  of  this 
figure  with  Figure  4  for  the  uniform  case  shows  that,  as 
expected,  tapered  illumination  results  in  power  density 
minima  instead  of  zeros.  Also,  the  envelope  peak  is  not 


constant  in  the  case  of  tapered  distributions  without 


pedestals . 


This  type  of  aperture  i 1 luninat ion  is  illustrated  in 
Figure  7  and  is  an  improvement  over  each  of  the  two 
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Figure  7.  Aperture  Illumination  with  Symmetric  Taper 
on  a  Uniform  Pedestal  (Sciambi,  1965:79) 


previously  discussed  distributions.  The  taper  accounts  for 
increased  attenuation  of  components  reflected  near  the 
paraboloid’s  edge,  and  the  pedestal  accounts  for  imperfect 
illumination  of  the  reflector  surface  by  the  feed  horn.  On 
author  has  said: 

Of  the  various  aperture  distributions  that  have 
been  studied,  one  in  particular  seems  to  satisfy 
most  needs  and  covers  typical  behavior  of  actual 
dish  antennas.  This  is  the  (1-P^  )**  on  a  ped¬ 
estal.  The  (l-p2)K  function  has  the  important 
advantage  that  both  pattern  and  directivity  can  be 
obtained  in  closed  form.  Inclusion  of  the  pedestal 
adds  versatility  and  allows  closer  representation 
of  actual  distributions.  (Hansen,  1964:64) 

The  generic  form  of  this  illumination  can  be  given  as 

G(x)  =  A  +  (l-A)(l-x*  )"  (11 

where 

0(x)  =  normalized  aperture  distribution 

A  =  amplitude  of  field  at  aperture  edge 
X  =  normalized  aperture  radial  variable 
N  =  exponent  controlling  rate  of  amplitude  taper 


The  distribution  has  two  independent  parameters,  A  and  N, 
and  is  normalized  to  one  at  the  aperture  center  where  x  =  0 
and  tapers  down  to  the  non-zero  value  of  A  at  the  aperture 
edge  where  x  =  1.  The  far  field  pattern  resulting  from  an 
illumination  of  this  form  can  be  obtained  by  adding  the 
patterns  produced  by  the  two  component  parts,  i.e.  the  taper 
and  the  pedestal  (Sciambi,  1965:79). 

Hansen  has  developed  a  distribution  which  is  different  in 
form  than  Eq  (11)  but  which  also  results  in  a  symmetric 
amplitude  taper  on  a  uniform  pedestal  (Hansen,  1975:184; 
Hansen,  1976b:50-52;  Hansen,  1 976a : 4 7 7-480  )  .  He  was 
dissatisfied  with  the  illumination  in  the  form  of  Eq  (11) 
because  it  "has  two  parameters,  and  there  is  no  simple  way 
of  finding  the  values  of  these  two  parameters  that  optimize 
one  quality,  e.g.  efficiency,  for  a  given  second  quality, 
e.g.  sidelobe  ratio"  (Hansen,  1975:184).  He  was  interested 
in  designing  an  illumination  model  that  would  permit 
circular  aperture  tradeoff  studies  based  on  a  single 
parameter.  The  result  of  his  work  was  the  one-parameter 
circular  aperture  distribution  of  the  form 

lo  (TH(  l-x2  )  ] 

G(x)  = -  (  12  ) 

lo  [TH] 

where 

G(x)  =  normalized  aperture  distribution 

lo  =  modified  Bessel  function,  zero  order 
=  normalized  aperture  radial  variable 
=  single  parameter  controlling  illumination 


X 

H 


This  distribution  is  normalized  to  one  at  the  aperture 
center  and  tapers  off  to  the  non-zero  value  of  l/Io(irH)  at 
the  aperture  edge.  The  H  parameter  is  related  to  the 
sidelobe  ratio  by 

2Ii  (irH  ) 

sir  =  17.57  dB  +  201og  - 

where 

sir  =  sidelobe  ratio 
Ii  =  modified  Bessel  function,  first  order 

It  can  be  seen  that  the  H  parameter  is  uniquely  related  to 
the  sidelobe  ratio;  once  it  is  determined  the  aperture 
distribution  is  known.  Beamwidth  and  aperture  efficiency 
can  also  be  determined  through  knowledge  of  the  H  parameter. 
Table  I  shows  the  theoretical  one-parameter  distribution 
values  as  determined  by  Hansen. 

This  one-parameter  aperture  distribution  has  been  used 
by  Hansen  (Hansen,  1976b:50-52)  to  calculate  Fresnel  region 
power  density.  He  concluded  that  the  resulting  curves 
"should  more  closely  represent  actual  parabolic  reflector 
and  circular  array  antennas"  than  those  resulting  from 
assuming  a  symmetric  taper  without  a  pedestal  and  that  "they 
are  appropriate  for  use  in  calculations  of  near  field 
radiation  hazards"  (Hansen,  1976b;52).  These  curves  are 
reproduced  in  Figures  8  and  9.  Note  that  in  comparison  to 
Figure  6,  where  the  aperture  distribution  tapered  to  zero  at 
the  edge,  inclusion  of  an  amplitude  pedestal  results  in 
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Figure  8.  Power  Density  vs  Range  for  the  One-parameter 
Aperture  Distribution  and  20/30  dB  Sidelobes 
(Hansen,  1976b:50) 


Figure  9.  Power  Density  vs  Range  for  the  One-parameter 
Aperture  Distribution  and  25/35  dB  Sidelobes 
(Hansen,  1976b:50) 
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IV. 


Design  and  Structure  of  POWERDEN 

Overview 

This  chapter  presents  the  details  of  how  POWERDEN 
(Appendix  A)  works  in  both  its  axial  and  off-axis  mode. 

Basic  relationships  for  determining  the  power  density  are 
discussed  along  with  the  Romberg  numerical  integration 
technique  which  is  used.  Limitations  on  aperture  size, 
axial  range,  and  off-axis  range  are  presented  followed  by  a 
discussion  of  the  program’s  input  requirements  and  output 
information  format.  The  chapter  concludes  with  a  discussion 
of  some  of  the  practical  aspects  of  using  POWERDEN. 

Basic  Relat ionships 

Power  Density .  If  an  observation  point  is  in  the  far 
field  or  Fresnel  region,  the  reactive  power  is  negligible, 
and  the  power  density  is  proportional  to  the  square  of  the 
field’s  magnitude  at  that  point  (Silver,  1949:113)  or 

S  ~  jEp|  2  (  14  ) 

whe  re 

S  =  power  density 
Ep  =  field  at  observation  point 

POWERDEN  calculates  the  power  density  according  to  Eq  (14) 
at  discreet  points  within  the  Fresnel  region  using  the  one- 
parameter  circular  aperture  distribution.  Since  this 
illumination  model  is  normalized  to  one,  the  actual  power 


density  at  these  points  is  obtained  by,  in  turn,  normalizing 
the  calculated  values  with  respect  to  a  known  quantity,  i.e. 
the  power  density  on  axis  at  a  range  of  20^ /X .  This 
distance  marks  the  beginning  of  the  far  field  so  its  value 
is  readily  determined  by  Eq  (1). 

Field  Strength ■  The  field’s  relative  magnitude  at  the 
observation  point  is  given  by 

I  I  1  +  cos ( 

Ep  ~  -  /  G{x)  Jo(ux)  exp -  x  dx  (15) 

'  A  •'o  8A 


where 


A  =  R/(2D2/X),  relative  range 
G(x)  =  lo [tH ( 1 -x^  )  ] ,  relative  aperture  illumination 
lo  =  modified  Bessel  function,  zero  order 
H  =  single  parameter  developed  by  Hansen 
Jo  =  Bessel  function,  zero  order 
u  =  kasin( 0 ) 

k  =  2ir/X ,  propagation  constant 
a  =  aperture  radius 

X  =  p/a,  normalized  aperture  radial  variable 


Eq  (15)  results  from  Eq  (6)  after  performing  the  integration 
with  respect  to  0'  and  making  the  substitutions  for  the 
normalized  aperture  radius  and  range  (Appendix  D).  The 
value  of  this  expression  is  obtained  by  numerical 
integration  and  then  squared  to  result  in  the  pre-normal ized 
value  of  power  density  as  discussed  in  the  last  section. 

H  Parameter .  The  H  parameter  is  related  to  the 
sidelobe  ratio  by  Eq  (13)  which  is  reproduced  here  in  a 
different  form: 


sir  -  17.57dB  -  201og 


[211  (TTH)! 


It  is  solved  for  by  a  subroutine  which  utilizes  the  Newton- 
Raphson  method  (Chapra  and  Canale,  1985:137-143)  for  finding 
the  root  of  Eq  (16).  An  initial  estimate  based  on  Hansen’s 
findings  (Table  I)  is  generated  internally  within  the 
program,  and  the  routine  converges  quadrat ically  for 
sidelobe  ratios  up  to  95  dB. 

Numerical  Integration 

The  real  and  imaginary  parts  of  Eq  (15)  are  integrated 
numerically  using  the  Romberg  adaptive  trapezoidal  technique 
(Chapra  and  Canale,  1985:416-424).  An  adaptive  method  was 
chosen  so  that  the  minimum  number  of  divisions  would  be  used 
at  each  observation  point.  The  integral  converges  more 
rapidly  for  points  near  the  axis  and  far  from  the  antenna; 
use  of  a  method  such  as  Simpson’s,  where  a  fixed  number  of 
divisions  is  used,  would  result  in  wasted  time  during 
evaluation  at  these  points.  The  number  of  divisions  is 
continually  doubled  up  to  a  maximum  of  4096  (2^^  )  or  until 
an  error  limit  set  by  the  operator  is  reached.  All  plots  in 
this  report  that  are  based  on  data  generated  by  POWERDEN 
were  arrived  at  using  the  default  error  limit  of  0.1% 


Limitations 

Aperture  Diameter .  The  method  developed  in  this 
report  is  valid  only  for  electrically  large  antennas.  The 
scalar  Kirchoff  diffraction  integral  (Eq  (3))  is  based  on 
this  assumption  as  is  the  design  of  paraboloids  by  the  laws 
of  geometric  optics.  POWERDEN  uses  a  value  of  ten 


wavelengths  for  the  minimum  aperture  diameter.  Most 
practical  paraboloids  are  much  larger  than  this  at  typical 
microwave  frequencies.  Below  this  limit  the  high  frequency 
assumptions  may  not  be  valid,  and,  at  close  range,  amplitude 
considerations  due  to  proximity  to  the  reactive  field  will 
dominate  phase  considerations  when  trying  to  establish  a 
minimum  axial  range  ( Eq  (7)). 

Axial  Range .  The  Fresnel  approximations  which  Eq  (15) 
is  based  on  are  not  valid  for  all  of  the  Fresnel  region. 

Just  as  the  far  field  of  a  large  paraboloid  begins  farther 
away  from  the  aperture  plane  than  that  of  a  smaller  one, 
the  minimum  range  at  which  the  Fresnel  approximations  are 
accurate  is  greater  for  a  large  antenna.  As  mentioned  in 
Chapter  2,  the  criterion  for  minimum  range  is 

R  >  0.5D(D/A)1/3  (7  ) 

In  the  axial  mode  of  operation  POWERDEN  calculates  power 
density  at  discreet  points  starting  at  the  beginning  of  the 
far  field  ( A  =  1 )  and  proceeding  toward  the  paraboloid  along 
the  boresight  until  the  minimum  distance  ( Eq  (7))  is  reached 
or  until  a  relative  range  of  A  =  0.01  is  reached,  whichever 
is  greater.  The  latter  condition  will  prevail  for  antennas 
wider  than  125x.  As  is  evident  from  an  examination  of 
Figures  8  and  9,  at  distances  less  than  A  -  0.01,  power 
density  is  varying  so  rapidly  with  respect  to  distance  that, 
in  terms  of  establishing  whether  the  power  density  is  at  a 
hazardous  level,  there  is  little  need  to  go  farther. 
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Off-axis  Range .  The  FSA  approximation  discussed  in 
Chapter  2  is  valid  for  angles  up  to  about  30°  off  the  axis 
of  symmetry  (Hansen,  1964:27,30,40).  This  range  covers  many 
sidelobes  for  a  large  antenna  near  the  far  field  and  covers 
the  first  four  sidelobes  of  a  ten  wavelength  antenna  with 
uniform  illumination  at  a  distance  of  2D^/\.  The  criterion 
used  by  POWERDEN  for  limiting  off-axis  range  was  given  in 
Chapter  2  as 


IT  £1  ^ 

-  sin2 (6)  <  0.08  (9) 

\z 

A  value  of  0.08  in  this  equation  allows  coverage  of  up  to 
26.8°  at  the  far  field  limit  of  2D^/\.  It  has  been  shown 
( Rahmat-Sami i  et  al . ,  1981:162)  that  the  FSA  approximations 
give  extremely  accurate  results  within  this  angular  range 
even  for  a  paraboloid  as  small  as  ten  wavelengths  in 
diameter . 

Program  Input 

Axial  Mode .  The  axial  mode  is  used  to  determine  if  a 
hazardous  power  density  level  exists  along  the  antenna 
boresight  and,  if  so,  what  the  maximum  value  is  and  at  what 
range  it  occurs.  In  this  mode  of  operation  POWERDEN 
requires  the  following  user  input: 

1.  Sidelobe  Ratio:  This  is  the  ratio  of  the  main 


beam  power  density  to  that  of  the  first  sidelobe. 


It  is 


given  in  decibels  and  is  limited  to  a  minimum  of  17.57  dB 
for  a  uniform  illumination  and  a  maximum  of  95  dB.  This 
maximum  value  is  much  larger  than  values  associated  with 
typical  paraboloidal  antennas.  The  sidelobe  ratio  is  used 
by  POWERDEN  to  calculate  the  H  parameter  so  that  the 
aperture  distribution  is  fixed  ( Eq  (12)). 

2.  Operating  Frequency:  Input  in  GHz. 

3.  Antenna  Diameter:  This  is  limited  to  a  minimum 
value  of  ten  wavelengths  and  is  input  in  units  of  meters. 

4.  Gain  and  Transmitted  Power:  These  values  are 
used  to  calculate  the  power  density  ( Eq  (1))  at  a  range  of 
2D2 /X .  This  figure  is  needed  for  determination  of  absolute 
power  densities  since  values  calculated  internal  to  POWERDEN 
are  normalized  with  respect  to  those  at  this  distance.  Gain 
is  input  in  decibels  and  power  in  kilowatts. 

5.  Maximum  Numeric  Integration  Error:  The  Romberg 
integration  technique  utilized  by  POWERDEN  compares  the  most 
recent  value  of  the  integral  with  that  computed  using  half 
as  many  divisions  during  the  previous  integration.  The 
number  of  divisions  continues  to  double  up  to  a  maximum  of 
4096  (2*2)  or  until  the  two  most  recently  calculated  values 
differ  by  no  more  than  this  error  limit.  The  default  error 
limit  is  set  at  0.1%. 

6.  Output  to  Line  Printer  or  Disc  File:  Data 


output  to  a  disc  file  might  be  desirable  either  if  the  line 
printer  is  being  used  or  if  it  cannot  be  allowed  to  stand 
idle  waiting  for  POWERDEN  to  complete  its  calculations.  If 


disc  file  output  is  selected,  POWERDEN  prompts  the  operator 
for  a  filename.  This  filename  must  be  a  valid  DOS  path. 

Of f-axi s  Mode .  The  off-axis  mode  is  used  for 
calculating  power  density  at  points  along  a  line 
perpendicular  to  the  axis  of  symmetry  starting  at  some  axial 
range  greater  than  the  limit  imposed  by  Eq  (7).  The  goal  is 
the  determination  of  a  hazard  radius  outside  of  which  the 
power  density  is  below  an  acceptable  level.  In  this  mode  of 
operation  POWERDEN  requires  all  of  the  same  inputs  just 
given  for  the  axial  mode  plus  the  following: 

1.  Off-axis  Starting  Point;  This  will  probably 
be  determined  from  information  obtained  using  POWERDEN  in 
the  axial  mode  as  well  as  the  potentially  hazardous 
situation  which  is  being  analysed.  The  input  is  given  in 
units  of  meters. 

2.  Off-axis  Stopping  Point:  The  optimum  value  of 
this  input  depends  strictly  on  the  individual  antenna  being 
analysed  and  the  experience  of  the  operator.  The  stopping 
point  will  be  a  value  which  can  be  assumed  by  the  expression 

u  =  kasin( $)  (17) 

where 

k  =  2ir/A,  propagation  constant 
a  =  antenna  radius 
6  =  off-axis  angle 


Construction  of  off-axis  radiation  patterns  in  terms  of  u 
suppresses  the  effect  of  diameter  size  so  that  the  patterns 
are  dependant  only  on  aperture  illumination  and  range.  The 


default  value  is  16,  and  this  will  cover  about  four 
sidelobes  at  the  beginning  of  the  far  field  of  a  uniformly 


illuminated  aperture.  If  power  density  at  the  off-axis 
starting  point  is  only  slightly  above  the  applicable 
standard,  it  probably  would  not  be  necessary  to  go  out  to  an 
angle  where  u  =  16.  On  the  other  hand,  an  extremely  high 
power  density  at  the  starting  point  might  require  a  larger 
angular  coverage  so  that  u  should  be  greater  than  16. 

Program  Output 

Axial  Mode .  Using  POWERDEN  in  the  axial  mode  results 
in  the  following  data  (Appendix  B)  being  output  to  a  line 
printer  or  disc  file: 

1.  Relative  Range:  POWERDEN  starts  at  a  relative 
range  of  ^  =  1  and  calculates  the  real  and  imaginary  parts 
of  Eq  (15)  for  196  points  along  the  axis  until  A  =  0.01. 

This  assumes  the  range  criterion  given  by  Eq  (7)  is  not 
violated.  These  points  are  not  evenly  distributed  between 
the  two  end  points;  smaller  steps  are  used  after  A  =  0.1  and 
even  smaller  ones  after  A  =  0.02.  This  greater  resolution 
insures  that  the  increasingly  rapid  variation  in  power 
density  at  closer  ranges  (Figures  8  and  9)  is  in  fact 
reflected  by  the  output  data. 

2.  Relative  Power  Density:  This  is  the  ratio 
expressed  in  decibels  of  power  density  at  the  particular 
point  on  axis  to  power  density  at  the  beginning  of  the  far 
field  where  A  =  1.  These  relative  power  density  figures 


make  the  detection  of  maximums  and  minimums  easier  and, 


in  conjunction  with  the  corresponding  figures  for  relative 

A 

range,  make  possible  the  production  of  plots  such  as  Figures  • 

8  and  9.  If  either  the  real  or  imaginary  part  of  the 
integral  in  Eq  (15)  fails  to  converge  at  a  point,  then  no 
relative  power  density  information  is  given  for  that  point. 

3.  Absolute  Range:  This  is  the  actual  distance  in 
meters  between  the  origin  of  the  aperture  plane  and  the 
point  at  which  power  density  has  been  calculated. 

4.  Absolute  Power  Density:  This  is  the  actual 
power  density  in  milliwatts  per  square  centimeter  at  the 
particular  axial  point.  As  with  the  relative  power  density, 
no  information  is  given  if  the  integral  did  not  converge  at  a 
particular  point. 

Off-axis  Mode .  Use  of  POWERDEN  in  the  off-axis  mode 
results  (Appendix  C)  in  the  output  of  absolute  power  density 
information  in  the  same  manner  as  with  the  axial  mode.  The 
other  output  is  the  following: 

1.  Value  of  kasin(^}:  POWERDEN  proceeds  off-axis 
in  increments  of  u  =  0.1  and  continues  at  this  same  rate 
until  the  limit  set  by  the  operator  is  reached. 

2.  Relative  Power  Density:  This  information 
differs  from  that  given  in  the  axial  mode  in  that  the  ratio 
is  with  respect  to  power  density  at  the  off-axis  starting 
point  rather  than  the  beginning  of  the  far  field.  This 
information  along  with  the  corresponding  values  of  u  makes 
possible  the  production  of  off-axis  radiation  plots  of  the 
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kind  introduced  in  the  next  chapter.  If  either  part  of  the 
integral  fails  to  converge  for  the  off-axis  starting  point, 
relative  power  density  cannot  be  given  for  any  of  the 
subsequent  points  since  all  of  them  are  normalized  with 
respect  to  the  value  at  this  starting  point.  Absolute  power 
density  data,  however,  will  still  be  available. 

3.  Distance  Off-Axis:  This  is  the  distance  in 
meters  from  the  field  point  to  the  antenna  boresight  along  a 
line  perpendicular  to  the  boresight.  This  information 
allows  determination  of  the  radiation  hazard  radius. 

Practical  Aspects 

On-screen  Information .  While  POWERDEN  is  running, 
information  about  the  calculations  at  each  data  point  is 
displayed.  A  typical  set  of  information  for  an  axial 
observation  point  might  be 

2  1.688104E-02 

4  -3.990152E-03 

8  2.073011E-02 

16  2.074877E-02 

Re  0.8400 

The  first  four  lines  show  the  number  of  divisions  used  by 
the  Romberg  integration  routine  and  the  corresponding 
integration  results.  Once  the  error  limit  set  by  the 
operator  has  been  reached,  the  last  line  is  displayed  and 
shows  that  this  calculation  is  for  the  real  part  of  the 
integral  at  a  relative  range  of  A  =  0.8400.  The  imaginary 
part  of  the  integral  is  computed  after  the  real  part  has 


been  determined  for  all  points  and,  in  that  case,  the 


letters  "Im"  are  displayed  instead  of  "Re".  An  on-screen 
information  set  is  the  same  for  the  off-axis  mode  except  the 
last  line  shows  the  value  of  u  instead  of  A. 

Running  Time .  Like  any  computer  program,  the  time 
required  by  POWERDEN  to  complete  a  set  of  calculations 
depends  to  a  large  degree  on  the  characteristics  of  the 
computer,  i.e.  clock  speed,  etc.  The  computer  used  to 
produce  the  data  used  for  this  report  was  a  government  owned 
Z-248  (AT  compatible)  equipped  with  a  numerical  math 
coprocessor  (80287  chip).  With  this  configuration,  typical 
run  times  for  both  axial  and  off-axis  modes  were  in  the  10- 
20  minute  range  depending  in  the  antenna  parameters  and 
numerical  integration  error  limit.  The  math  coprocessor  is 
highly  recommended  since  its  absence  can  result  in  a  tenfold 
increase  in  run  time. 


V. 


Results 

Overview 

This  chapter  presents  in  graphical  form  the 
applying  POWERDEN  to  the  analysis  of  a  uniform 
illumination  (H  =  0)  and  to  four  others  in  the 
one-parameter  distribution  ( Eq  (12))  with  assumed  sidelobe 
ratios  of  20,  25,  30,  and  35  dB.  Axial  results  are  compared 
directly  with  those  achieved  by  Hansen.  Off-axis  results 
for  the  one-parameter  distribution  have  not  been  published 
so  direct  comparison  is  not  possible  except  in  the  far  field 
case  of  the  uniform  distribution  since  it  is  available  in 
most  books  dealing  with  paraboloidal  analysis.  Therefore, 
off-axis  results  for  the  four  other  tapered  distributions 
will  be  presented  and  analysed  in  relation  to  what  is 
expected  theoretically.  The  goal  will  be  to  establish  that 
POWERDEN  yields  reasonable  data  and  can  be  a  useful  tool  for 
determining  hazardous  radiation  zones  in  the  Fresnel  region. 

Axial  Data 

Figures  10  and  11  show  the  results  of  using  POWERDEN  in 
the  axial  mode  for  the  analysis  of  the  four  tapered 
distributions.  Direct  comparison  with  Hansen’s  published 
results  (Figures  8  and  9)  shows  that  POWERDEN  achieves 
nearly  identical  results.  The  less  tapered  distributions 
with  the  highest  sidelobes  undergo  the  most  significant  near 
field  gain  reduction.  Power  density  fluctuations  begin 


results  of 
aperture 
form  of  the 


Figure  10.  POWERDEN  Axial  Radiation  Pattern  for  Sidelobe 
Ratios  of  20/30  dB 


Figure  11.  POWERDEN  Axial  Radiation  Pattern  for  Sidelobe 
Ratios  of  25/35  dB 


around  a  relative  range  of  A=  0.1  as  the  aperture  begins 
to  encompass  an  increasing  number  of  Fresnel  zones.  The 
magnitude  of  the  fluctuations  decreases  as  the  illumination 
becomes  more  highly  tapered;  this  is  due  to  the  decreasing 
field  strength  of  the  outer  Fresnel  zones  in  a  tapered 
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illumination  compared  to  the  inner  ones.  In  contrast  to 
Figure  6  where  the  aperture  distribution  went  to  zero  at 
the  paraboloid’s  edge,  i.e.  no  pedestal,  here  the 
oscillation  peaks  are  nearly  constant  for  a  given  sidelobe 
ratio.  (Hansen,  1976b:51) 


Off-Axis  Data 

Far  Field .  Radiation  patterns  resulting  from  use  of 
POWERDEN  in  the  off-axis  mode  at  a  single  range  of  2D^ /a  are 
presented  in  this  section  so  that  the  effect  of  increasing 
illumination  taper  can  be  seen  without  being  obscured  by 
near  field  phenomena.  Increasing  taper  should,  in  general, 
result  in  wider  beamwidth,  lower  gain,  and  lower  sidelobe 
levels  (Silver,  1949:179;  Jull,  1981:22). 

Figure  12  is  the  theoretical  far  field  radiation  pattern 
of  a  uniform  distribution  while  Figure  13  was  produced  from 
POWERDEN  output  data.  It  can  be  seen  that  the  two  figures 
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Figure  14. 


POWERDEN  Far  Field  Radiation  Pattern 
(sir  =  20  dB) 


Relative 
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Figure  16  displays  a  sidelobe  level  cf  about  29  dB  instead 
of  the  30  dB  which  was  input.  Beamwidths  are  also  slightly 
larger  than  those  given  in  Table  I.  The  deviations  from 
theoretical  expectations  are  quantified  in  Table  II.  The 
POWERDEN  sidelobe  ratios  are  those  resulting  from  the 


Table  II 

Comparison  of  Theoretical  One-Parameter  Distribution  Values 


With  Those  Resulting  From  POWERDEN 


I  nput 

POWERDEN 

Beamwidth 

S  idelobe 

Sidelobe 

From 

POWERDEN 

Ratio 

Ratio 

Table  I 

Beamwidth 

(dB) 

(dB) 

{ radians  ) 

( radians  ) 

17.57 

17.31 

1.0290  A/D 

1.0295  a/d 

20 

19.63 

1.0787 

1.0800 

25 

24 . 30 

1.1734 

1,1758 

30 

28 . 78 

1.2607 

1 . 2635 

35 

33.00 

1 . 3403 

1 . 3437 

resolution  provided  by  incrementing  the  angular  range  in 
steps  of  u  =  0.1,  and  the  POWERDEN  beamwidths  were  computed 
by  linear  interpolation  between  the  two  values  of  u  which 
fell  on  either  side  of  -3.01  dB.  Although  the  sidelobe 
levels  are  elevated,  the  practical  effect  will  be  the 
establishment  of  a  wider  radiation  hazard  radius  than  is 
absolutely  necessary.  This  same  observation  can  be  made 
regarding  the  slightly  wider  beamwidths;  the  theoretical 
half-power  point  is  closer  to  the  antenna  boresight  so  the 
radiation  hazard  radius  will  be  slightly  larger  than  called 


for  in  the  ideal  case. 


Fresnel  Region .  The  quadratic  phase  term  in  Eq  (15) 

was  included  so  that  the  phase  criterion  of  t/8  used  in  the 

far  field  could  be  maintained  in  the  Fresnel  region.  The 

effect  of  finite  range  has  been  described  as  follows: 

In  the  region  closer  than  20^ /X ,  called  the 
radiating  near  field  region,  the  pattern  varies 
with  the  distance  and,  due  to  the  quadratic  phase, 
perfect  interference  does  not  occur.  This  blurs 
the  nulls  and  raises  the  sidelobes.  At  closer 
distances  ...  the  main  beam  broadens  and  some 
close-in  sidelobes  may  be  swallowed  by  the  main 
beam.  (Hansen,  1964:29) 

In  addition,  the  quadratic  phase  has  a  greater  effect  on  the 
first  null  than  on  subsequent  ones  (Rusch  and  Potter, 
1970:99),  and  it  has  a  smaller  overall  effect  on  the  more 
highly  tapered  distributions  (Jull,  1981:32). 

Figures  18  and  19  illustrate  these  effects  for  the  case  of 
a  rectangular  aperture  with,  respectively,  uniform  and 
tapered  illumination.  The  scales  are  different  than  those 
used  for  other  figures  in  this  report,  but  they  do 
illustrate  the  main  effects  of  finite  range  (  /3  =  0  is  at  a 
relative  range  of  A  =  ®  while  0  =  ir/2  is  at  a  relative  range 
of  A  =  0.25).  The  sidelobes  are  elevated  for  A  =  0.25,  and 
each  figure  shows  a  broadened  main  beam  which,  in  Figure  19 
for  the  tapered  distribution,  has  completely  swallowed  the 
first  sidelobe.  In  addition,  the  first  null  in  each  far 
field  plot  is  blurred  more  than  the  second. 

Figures  20-24  show  the  off-axis  results  achieved  by 
POWERDEN  in  the  Fresnel  region  of  the  uniform  and  tapered 
distributions.  For  each  illumination  results  are  plotted  for 
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relative  ranges  of  A  =  1,  A  =  0.5,  and  A  =  0.25.  Note  that, 
as  predicted,  the  sidelobes  are  elevated  with  decreasing 
range  and  the  beamwidths  are  broadened.  At  A  =  0.25  the  30 
and  35  dB  sidelobes  are  nearly  totally  engulfed  by  the  main 
beam.  In  each  of  the  five  cases  the  first  null  is  affected 
more  significantly  than  is  the  second.  The  more  highly 
tapered  distributions  with  the  larger  sidelobe  ratios  suffer 
less  overall  distortion  from  the  quadratic  phase  effects  of 


the  closer  ranges. 
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VI .  Conclusion 


Overview 


The  results  presented  in  the  last  chapter  show  that  the 


data  produced  by  POWERDEN  in  the  axial  mode  corresponds 


almost  exactly  with  that  published  by  Hansen.  The  off-axis 


results  cannot  be  compared  directly  since  none  have  been 


published  for  the  one-parameter  circular  aperture 


distribution.  However,  the  data  produced  by  POWERDEN  in  the 


off-axis  mode  is  consistent  in  all  respects  with  that  which 


would  be  expected  theoretically. 


Recommendations  for  Further  Research 


SBSSi 


The  goal  of  this  thesis  was  the  development  of  a  better 
method  for  Air  Force  use  in  determining  the  Fresnel  region 


power  density  hazards  associated  with  paraboloidal  antennas. 
POWERDEN  accomplishes  that  goal;  it  corrects  the  three 
deficiencies  of  the  present  method  discussed  in  Chapter  1 
by:  (1)  incorporating  a  more  realistic  aperture  illumination 
model,  (2)  being  applicable  to  a  spectrum  of  specific 
illuminations  instead  of  the  most  satisfactory  of  four 
predetermined  distributions,  and  (3)  unambiguously 
estimating  the  specific  illumination  from  the  sidelobe  ratio 
instead  of  the  beamwidth.  The  program  is  meant  to  be  a 
useful  tool  for  use  by  those  already  familiar  with 
paraboloidal  radiation  characteristics  and  the  related 
problems  of  electromagnetic  compatibility. 


V.V.V.  V.  V. 


els 

print 

print 

print 

print 

print 

print 

print 

print 

print 

print 

print 

print 

print 

print 

print 

print 

print 

print 

print 

print 

print 

print 


Appendix  A:  POWERDEN  Source  Code 


print 

1lHHHt.*iHHHliif%1f$f************  ****************' 
*  ♦' 

*  POWERDEN  ♦' 


♦  Johnnie  E.  Mize,  Capt ,  USAF  ♦" 

*  Air  Force  Institute  of  Technology  ♦” 

♦  Class  GE-87D  *" 

*  ♦ " 

POWERDEN  is  designed  to  calculate  the  Fresnel" 
region  power  densities  associated  with  high  gain," 
center  fed,  paraboloidal  antennas.  It  assumes  a" 
diameter-to-wavelength  ratio  of  at  least  ten." 

The  illumination  is  modeled  as  a  symmetric  taper" 
on  a  uniform  pedestal  with  constant  phase  and" 
linear  polarization." 
print 


The  aperture  illumination  model  used  in  this  program 
is  the  one-parameter  circular  distribution  developed 
by  R.  C.  Hansen.  Details  have  been  published  in: 

1)  IEEE  Transactions  on  Antennas  and  Propagation, 

AP-24:  477-480  (July  1976) 

2)  Microwave  Journal,  19:  50-52  (February  1976) 

3)  Electronic  Letters,  11:  184  (April  1975) 

The  program  uses  Romberg  integration  to  solve  the  scalar 
Kirchoff  diffraction  integral  in  the  Fresnel  Region. 

A  glossary  of  the  main  variables  is  located  at  the  end 
of  this  source  code. 


input  "Press  <RETURN/ENTER>  to  proceed.  ",ansX 
els 

Begin: 

input  "Input  the  first  sidelobe  ratio  (17.57  -  95  dB ) :  ",slr 
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.w'  • 


*k '  •% 


if  sir  >=  17.56999  and  sir  <=  95  then 
print 

row  =  csrlin 
print  "Working  ..." 
end  if 

il#  =  (10^((8lr  -  1 7 . 5 7  ) /20  )  ) /2  ’  from  Hansen’s  formula 

pi#  =  4*atn( 1  ) 

tolerance  =  lE-14  ’  accuracy  of  internal  calculations 

if  sir  <  17.56999  then  ’  choose  initial  estimate  for  using 

print  "Too  small."  ’  Newton’s  method  to  solve  for  the 
print  ’  H  parameter 

goto  Begin 

elseif  sir  <  17.57001  then 
H#  =  0 
goto  Zero 

elseif  sir  <=  25  then 
est#  =  2.9 

elseif  sir  <=  50  then 
est#  =  7 

elseif  sir  <=  75  then 
est#  =  10.4 
elseif  sir  <=  95  then 
est#  =  13.1 
else 

print  "Too  large." 
print 

goto  Begin 
end  i  f 

’  calculate  the  H  parameter 
do  until  fnF#(est#)  <=  tolerance 

est#  =  est#-fnF# ( est# ) /fnG# { est# )  ’ Newton-Raphson  formula 
1  oop 

H#  =  est#/pi# 

Zero  : 

locate  row , 1 

print  "The  H  parameter  is:"; 
print  using  "  #.####";H# 

Recycles  : 


print 

input  "Input  the  operating  frequency  (GHz):  ",freq# 

if  freq#  <=  0  then  print  "Too  small.”  :  goto  Recycles 

lambda#  =  0.30/freq# 
k#  =  2*pi#/lambda# 
print 

print  "The  wavelength  is" ; 100*lambda# ;  "centimeters." 
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Recyclel  : 


print 

input  "Input  the  antenna  diameter  (meters):  ",diam# 

if  diam#<9 . 99*lambda#  th^n  print  "Too  small."  :goto  Recyclel 

RangeMin#  =  0 . 5*d iam#* ( d iam#/l ambda# ) * ( 1 /3 ) 
ffmin#  =  2*diam#*2/lambda# 
a#  =  diam#/2 
print 

print  "The  far  field  begins  at  a  range  of f fmin# ; "meters . " 

Recycle9 : 

print 

input  "Input  the  antenna  gain  and_ 

transmitted  power  (dB,kW):  ", gain, power 

if  gain  <=  0  or  power  <=  0  then_ 

print  "Too  small."  :  goto  Recycle9 

ffpwrden#  =  1 00*power» { 1 0* ( 0 . 1 ♦gain ) ) /_ 

( 4»pi#*f fmin# " 2 )  ’  mW/cm“2 

print 

print  "Power  density  at  the  beginning  of  the  far  field  is"; 
print  using  "  #.######""'"  ";ffpwrden#; 
print  "mW/cm*2." 

Recycle2 : 

print 

input"Do  you  need:  on-axisll)  or  off-axi8(2)  data?  ",ansX 

if  ansX  =  1  then 
exit  if 

elseif  ansX  =  2  then 
goto  OffAxisStart 
else 

goto  Recycle2 
end  i  f 

OnAxi sStart : 

print 

key  10,  "0.1"  +  chr$(13) 

print  "Input  max  error  for  numeric  integration"; 
input  "  (Default  =  O.IX,  Press  FIO):  ",ErrLimOn 

if  ErrLimOn  =  0  then  print  "Too  small."  :  goto  OnAxisStart 

key  10,  chr$(00) 

ErrLimOn  =  ErrLimOn/100 


print 

input  "Write  data  to:  line  printer{l)_ 
or  disc  file(f)  ?  ", answer! 

if  answer!  =  "1"  then 
exit  if 

elseif  answer!  =  "f"  then 
print 

input  "Filename:  ",B! 
else 

goto  RecyclelO 
end  i  f 

print 

input  "Press  <RETURN/ENTER>  to  proceed.  ",anB% 
print  :  print 

if  DimFlagSX  =  1  then  DimFlagSX  =  0  :  erase  RealOn# , ImOn# 

dim  RealOn# ( 0 : 1 95 ) , ImOn# ( 0 : 1 95 ) ’ to  store  integration  results 

DimFlagSX  =  1 

flagX  =  0  :  A!  =  "  Re" 

u  =  0  ’so  that  fnJ#(x#)  =  1 

FirstDecade:  ’  relative  range  1  to  0.1 

aX  =  100  :  bX  =  10  :  dX  =  -1  :  eX  =  10 
elementX  =  -1 

NextDecade :  ’  relative  range  0.1  to  0.02,  0.02  to  0.01 

eX  =  10*eX 

for  RX  =  aX  to  bX  step  dX 
R  =  RX/eX 
incr  elementX 

if  DimFlaglX  =  1  then  DimFlaglX  =  0  :  erase  Romberg# 

dim  Romberg# ( 1 : 1 2 , 1 : 1 2 )’ store  intermediate  Romberg  values 
DimFlaglX  =  1 

’  check  for  distance  from  aperture 
if  R*ffmin#  <  RangeMin#  then  goto  StartImOn 
q#  =  pi#/(8*R) 

’  compute  initial  Romberg  seed 

if  flagX  =  0  then 

Romberg# ( 1 , 1 ) =0 . 25* ( cos ( q# ) +2*f nMidSum# ( 2 ) ) 

else 

Romberg#  (l,l)=0.25*{8in(q#  )  ■•■2*f  nMidSum#  (  2  )  ) 
end  i f 


pi int  using  "#### 


"  :  nX : 


if  Romberg# ( 1 , 1 )  <  0  then 

print  using  Romberg#  (  1  ,  1  ) 

else 

print  using  "  #.######"^  Romberg# ( 1 , 1 ) 

end  i  f 

SubDivideOn:  ’  more  divisions  for  greater  accuracy 

nX  =  2*nX  :  jX  =  jX  +  1  :  kX  =  1 

’  exit  if  no  far  field  convergence 
if  nX  =  8192  and  elementX  =  0  then 
print 

print  "No  Convergence  In  The  Far  Field." 
goto  Recycle4 
end  i  f 

’  exit  loop  after  2*12  divisions 
if  nX  =  8192  and  flagX  =  0  then 
RealOn# ( elementX )  =  0 
goto  MarkO 

elseif  nX  =  8192  and  flagX  =  1  then 
ImOn# ( elementX )  =  0 
goto  MarkO 
end  i  f 

’  compute  subsequent  Romberg  seeds 

if  flagX  =  0  then 

Romberg# (jX, 1 )  =  (0. 5/nX)*(cos( q# )  +  2*fnMidSum#(nX)  ) 

else 

Romberg# (jX, 1 )  =  (0.5/nX)*(sin(q# )  +  2*f nMidSum# ( nX  )  ) 
end  i  f 

’  calculate  intermediate  Romberg  values 
for  qX  =  (jX-1)  to  1  step  -1 
kX  =  kX  +  1 
w#  =  4* ( kX-1  ) 

Romberg# ( qX , kX )  =  ( w#* Romberg# ( qX+ 1 , kX- 1  )_ 

-Romberg# ( qX , kX- 1 )  )/(w#-l  ) 

next  qX 

print  using  "####  " ; nX ; 

if  Romberg# ( 1 , kX  )  <  0  then 

print  using  "+#.######  Romberg# ( 1 , kX ) 

else 

print  using  "  #.######*** ~"; Romberg# ( 1 , kX ) 
end  i  f 

er#  =  1  -  Romberg# ( 1 , kX- 1 ) /Romberg# ( 1 , kX ) 

’  check  for  max  error 
if  abs(er#)  <  ErrLimOn  and  flagX  =  0  then 
RealOn# ( elementX )  =  Romberg# ( 1 , kX ) 
elseif  abs(er#)  <  ErrLimOn  and  flagX  =  1  then 
ImOn# ( elementX  )  =  Romberg# ( 1 , kX ) 
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goto  SubDivideOn 
end  i  f 

MarkO:  ’  display  info  about  current  data  point 

print  A$ ; 

print  using  ”  #.####";R; 

if  nX  =  8192  then 

print  "  No  Convergence  At  This  Error  Limit" 

end  i  f 

print  :  print  :  print 
next  R% 


if  eleme 

ntX 

=  90  then  aX  =  99  :  bX  = 

20  ;  dX  =  -1 

:  goto_ 

NextDecade 

if  eleme 

ntX 

=  170 

then  aX  =196  :  bX 

=  100  :  dX  = 

-4  :  goto_ 

NextDecade 

StartImO 

n  : 

if  flagX 

= 

1  then 

goto  CalculateOn 

flagX  = 

1  ; 

A$  =  " 

Im"  :  goto  FirstDecade 

CalculateOn 

’  calculat 

e  and  display 

axial  data 

if  answe 

r$ 

=  "f"  then  goto  DiscWrite 

On 

print 

print  "J 

ob 

complet 

e  ..." 

Ipr int 

Iprint  " 

Axial  Power 

Density  Data" 

Iprint  : 

Ip 

r  int 

Iprint 

Relative 

Absolute " 

Iprint  " 

Power 

Absolute 

Power" 

Iprint  " 

Re 

lat  i  ve 

Dens ity 

Range 

Dens ity " 

Iprint  " 

Range 

( decibels  ) 

( meters  ) 

( mW/cm" 2  )  " 

Iprint 

elementX 

= 

-1 

norm#  = 

(  Re 

alOn#(0 

)"2  +  IiDOn#(0)"2) 

’  normal  i  ze 

wrt  ffmin# 

aX  =  100 

bX  =  10 

:  dX  =  -1  :  eX  = 

10 

Marks : 

eX  =  10* 

eX 

for  RX  = 

aX 

to  bX 

step  dX 

R  =  RX/eX 


if  R^ffmin#  <  RangeMin#  then  goto  ParameterListOn 


incr  elementX 

s#  =  (  RealOn#  (  elementX  ) ''2  +  ImOn#  {  elementX  )  "2  )/(  R"2*norm#  ) 
Iprint  using  "  #.####";R; 

if  RealOn# ( elementX )  =  0  or  ImOn# ( elementX )  =  0  then 
Iprint  ”  ♦***»»*". 

else 

Iprint  using  "  #######.####"; 1 0*log 1 0 ( s#  )  ; 

end  if 

Iprint  using  "  #.######"'  ";R*ffmin#; 

if  Real  On# ( el ementX  )  =  0  or  ImOn# ( elementX )  =  0  then 
Iprint  "  ♦♦♦**»***♦♦»" 

else 

Iprint  using  "  s#* f fpwrden* 

end  i  f 
next  RX 

if  elementX  =  90  then  aX  =  99  : bX  =  20  : dX  =  -1  :  goto  Mark5 

if  elementX  =  170  then  aX  =  196  :bX=100  : dX  =  -4  :goto  MarkS 
ParameterListOn : 

Iprint 

Iprint  "  H  parameter  ="; 

Iprint  using  " 

Iprint  "  sidelobe  ratio  ="; 

Iprint  using  "  ##.###  ";slr; 

Iprint  "dB" 

Iprint  "  frequency  ="; 

Iprint  using  ”  #.#####^''"^  ";freq#; 

Iprint  "GHz" 

Iprint  "  diameter  ="; 

Iprint  using  "  I.#####''""  ";diam#; 

Iprint  "meters” 

Iprint  "  gain  =" ; 

Iprint  using  "  ##.####" ^ " "  ";gain; 

Iprint  "dB" 

Iprint  "  transmitted  power  =" ; 

Iprint  using  "  #.######"'"''  ";power; 

Iprint  "kW" 

Iprint  "  error  limit  ="; 

Iprint  using  "##.###"; 1 00*ErrLimOn ; 

Iprint  "X" 
goto  Recycles 

DiscWriteOn ; 

print 

print  "Job  complete;  Writing  data  to  file 
print  B$ 

open  B$  for  output  as  #1 


print  #1  , 

print  #1,"  Axial  Power  Density  Data" 

print  #1,  ;  print  #1, 

print  #1,"  Relative 

print  #1,"  Power  Absolute 

print  #1,"  Relative  Density  Range 

print  #1,"  Range  (decibels)  (meters) 

print  #1, 


Absolute" 

Power" 

Density" 

(mW/cm"2)" 


elementX  =  -1 

norm#  =  ( RealOn# ( 0 ) " 2  +  ImOn#(0)"2)  ’  normalize  wrt  ffmin# 

aX  =  100  :  bX  =  10  :  dX  =  -1  :  eX  =  10 

f . Marks : 

eX  =  10*eX 

for  RX  =  aX  to  bX  step  dX 
R  =  RX/eX 

if  R*ffmin#  <  RangeMin#  then  goto  f . ParameterListOn 
incr  elementX 

s#  =  (  RealOn#  (  elementX  )  ^  2  +  ImOn#  (  elementX  )'' 2  )/(  R' 2*norm#  ) 
print  #1,  using  "  #.####";R; 

if  RealOn# ( elementX )  =  0  or  ImOn# ( elementX )  =  0  then 
print  #1,  " 
else 

print  #1,  using  "  #######.####"; 1 0*log 1 0 ( s#  )  ; 

end  i  f 


print  #1 ,  using  " 


«.####«#  " ;R*ffmin# ; 


if  RealOn# ( elementX )  =  0  or  ImOn# ( elementX )  =  0  then 
print  #1,  " 
else 

print  #1,  using  "  #•######""''''";  s#* f  fpwrden# 

end  i  f 


next  RX 


if  elementX  =  90  then  aX  =  99  :bX=20  : dX  =  -1  :goto  f.MarkS 
if  elementX  =  170  then  aX=l 96 : bX= 1 00 : dX=-4 : goto  f.MarkS 
f . ParameterListOn : 


print  #1, 

print  #1,  "  H  parameter  ="; 

print  #1,  using  "  #.####" ;H# 

print  #1,  "  sidelobe  ratio  ="; 

print  #1,  using  "  ##.###  ";slr; 
print  #1,  "dB" 

print  #1,  "  frequency  ="; 
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using  "  #.##### . ;freq#; 

"GHz" 

"  diameter  ="; 

using  "  #.#####  ";diam#; 

"meters" 

"  gain  ="; 

using  "  ##.####“"""  ";gain; 

"dB" 

"  transmitted  power  =" ; 

using  "  #  .######"'^"''  ";power: 
"  kW" 

error  limit  =  ; 
using  "##.###";  1 00*ErrLiiDOn  ; 

tt  a/  n 


print  #1  , 
print  #  1  , 
print  #1  , 
print  #1  , 
print  #1, 
print  #1 , 
print  #1, 
print  #1, 
print  #1 , 
print  #1, 
print  # 1 , 
print  #1 , 
print  # 1 , 
print  #1 , 
close  #1 

Recycles  : 


print 

input  "Do  you  need  off-axis  data  U8ing_ 

the  current  parameters  (y/n)?  ",ans$ 

if  ans$  =  "y"  then 
goto  OffAxisStart 
elseif  ans$  =  "n"  then 
exit  if 
else 

goto  Recycles 
end  i  f 


Recycle4  : 


print 

input  "New  problem  (y/n)?  ",ans$ 

if  ans$  =  "y"  then 
els 

goto  Begin 

elseif  ans$  =  "n"  then 
exit  if 
else 

goto  Recycle4 
end  if 
end 

Of  f Ax i sStart : 


print 

input  "Start  off-axis  calculations  at  what_ 
axial  range  (meters)?  ",z# 

if  z#  <  RangeMin#  then 
print  "Too  close." 
goto  OffAxisStart 


elseif  z#  >  ffmin#  +  0.001  then 

print  "That’s  in  the  far  field." 
goto  OffAxisStart 
end  i  f 

*  compute  max  angle  off-axis 

SinSquareThetaMax#  =  { 0 . 08*lambda#*z# ) / ( pi#*a# ^ 2 ) 

ThetaMax#  =  atn( sqr ( SinSquareThetaMax#/_ 

(1  -  SinSquareThetaMax#))) 

Recycle? : 

key  10,  "16"  +  chr$(13) 
print 

print  "Input  a  value  (Default  =  16,  Press  FIO)  for_ 
the  range  of  kasin("; 
print  chr$ ( 233  )  ; 
input  " ) :  " , f 

if  abs(f)  >  k#*a#*sin( ThetaMax# )  then 
print  "Too  large;  max  ="; 
print  k#*a#»sin( ThetaMax# ) 
goto  Recycle? 
else 

f%  =  fix(10*f) 
key  10,  chr$(00) 


Recyclel 2 ; 
print 

key  10,  "0.1"  +  chr$(13) 

print  "Input  max  error  for  numeric  integration"; 
input  "  (Default  =  O.IX,  Press  FIO):  ",ErrLimOff 

if  ErrLimOff  =  0  then  print  "Too  small."  :  goto  Recyclel2 

key  10,  chr$(00) 

ErrLimOff  =  ErrL imOf f / 1 00 

Recyclel 1 : 

print 

input  "Write  data  to;  line  printer(l)_ 
or  disc  file(f)  ?  ".answer# 

if  answer#  =  "1"  then 
exit  if 

elseif  answer#  =  "f"  then 
print 

input  "Filename:  ",B# 
else 

goto  Recyclell 
end  i  f 


print 

input  "Press  <RETURN/ENTER>  to  proceed.  ",ans% 
print 

flagX  =  0  :  A$  =  "  Re" 

if  DiBiFlag2X  =  1  then  DimFlag2X  =  0  :  erase  RealOf  f # ,  ImOf  f  # 

’to  store  integration  results 
dim  RealOff#(0:fX) ,ImOff#(0:fX) 

DimFlag2X  =  1 

Markl  : 

NormFlagX  =  0 
print 

row  =  csrlin  -  1 
print  "Working  ..." 

Mark2 : 

eiementX  =  -1 

for  uX  =  0  to  fX  ’  range  of  ka8in(theta) 

u  =  uX/10 
incr  elementX 

if  DimFlaglX  =  1  then  DimFlaglX  =  0  :  erase  Romberg# 

dim  Romberg# ( 1 : 1 2 , 1 ; 1 2 )’ store  intermediate  Romberg  values 
DimFlaglX  =  1 

SinSquareTheta#  =  ( u/( k#*a# ) ) "2 

Theta#  =  atn( Bqr( SinSquareTheta#/ ( 1-SinSquareTheta# ))  ) 

if  NormFlagX  =  0  then  R#  =  1_ 

else  R#  =  z#/ ( f fmin#*cos ( Theta# )  ) 

q#  =  pi#/(8*R#) 

’  compute  initial  Romberg  seed 

if  flagX  =  0  then 

Romberg# ( 1 , 1  )  =  0 . 25* ( cos ( q# ) *f nJ# ( u )  +  2* f nMidSum# ( 2 )  ) 
else 

Romberg# ( 1 , 1  )  =  0 . 25* ( s in( q# ) * f nJ# ( u )  +  2* f nMidSum# ( 2 ) ) 
end  i  f 

nX  =  2  :  jX  =  1 

if  NormFlagX  =  1  and  elementX  =  0  then  locate  row,l 

if  NormFlagX  =  1  then 

print  using  "####  ";nX; 

if  Romberg# ( 1 , 1 )  <  0  then 

print  using  "+#.######  Romberg# ( 1 , 1 ) 

else 

print  using  "  #.######*  Romberg# ( 1 , 1 ) 

end  i  f 
end  i  f 
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SubDivideOf  f : 


more  divisions  for  greater  accuracy 


nX  =  2*nX  :  jX  =  jX  +  1  :  kX  =  1 

’  exit  if  no  far  field  convergence 

if  nX  =  8192  and  NormFlagX  =  0  then 
locate  row , 1 

print  "No  Convergence  In  The  Far  Field." 
goto  Recycle6 
end  i  f 

if  nX  =  8192  and  flagX  =  0  then 

RealOf  f #  ( elementX )  =  0  ’exit  loop  after  2''12  divisions 
goto  Mark3 

elseif  nX  =  8192  and  flagX  =  1  then 
ImOf f# ( elementX )  =  0 
goto  Marks 
end  if 

’  compute  subsequent  Romberg  seeds 

if  flagX  =  0  then 

Romberg# ( jX , 1  )  =  (0.5/nX)*_ 

(coa(q#)»fnJ#(u ) +2*fnMidSum# ( nX ) ) 

else 

Romberg# ( jX , 1 )  =  (0.5/nX)*_ 

(sin(q#)*fnJ#(u )  +  2 *f nMidSum# ( nX )  ) 

end  i  f 

’  calculate  intermediate  Romberg  values 
for  qX  =  (jX-1)  to  1  step  -1 
kX  =  kX  +  1 
w#  =  4" (kX-1  ) 

Romberg# ( qX , kX ) =( w#* Romberg# ( qX+ 1 , kX- 1 )_ 

-Romberg# ( qX, kX-1 )  )/(w#-l  ) 

next  qX 

if  NormFlagX  =  1  then 

print  using  "####  '';nX; 

if  Romberg# ( 1 , kX )  <  0  then 

print  using  "+#.###♦##  ^ Romberg# ( 1 , kX ) 

else 

print  using  "  #.######  Romberg# ( 1 , kX ) 

end  i  f 
end  i  f 

er#  =  1  -  Romberg# ( 1 , kX-1 ) /Romberg# ( 1 , kX ) 

’  check  for  max  error 

if  abs(er#)  <  ErrLimOff  and  flagX  =  0  then 
RealOf f# ( elementX )  =  Romberg# ( 1 , kX ) 
elseif  abs(er#)  <  ErrLimOff  and  flagX  =  1  then 
ImOff# ( elementX )  =  Romberg# ( 1 , kX ) 
else 

goto  SubDivideOff 
end  i  f 


V 


Marks : 


display  info  about  current  data  point 


if  NormFlagX  =  1  then 
print  A$ : 

print  using  "  #####.#" ;u; 

if  nX  =  8192  then 

print  "  No  Convergence  At  This  Error  Limit 

end  if 

print  :  print  :  print 
end  i  f 


’  store  normalization  data 
if  NormFlagX  =  0  and  flagX  =  0  then 
HoldRealNorm#  =  Real0ff#{0) 

NormFlagX  =  1 
goto  Mark2 

elseif  NormFlagX  =  0  and  flagX  =  1  then 
HoldImNorm#  =  ImOff#(0) 

NormFlagX  =  1 
goto  Mark2 
end  i f 


next  uX 


StartImOf  f ; 


if  flagX  =  1  then  goto  CalculateOff 
flagX  =1  :  A$  =  ’’  Im"  :  goto  Markl 

CalculateOff:  ’  calculate  and  display  off-axis  data 

if  answers  =  "f"  then  goto  DiscWriteOff 


print 

print 

Iprint 

Ipr int 

Iprint 

Iprint 

Iprint 

Iprint 

Iprint 

Iprint 

Iprint 

Iprint 


"Job  complete 


:  Iprint 

tt 


"  kasin( " ; 
chr$(233  )  ; 


Off-Axis  Power  Density  Data" 


Relative 

Power 

Density 


( decibels ) 


Distance 
Of  f-Ax i s 

( meters ) 


Absolute" 

Power" 

Density" 


( mW/cm*2 ) 


elementX  =  -1  ’  normalize  wrt  ffmin# 

norm#  »  4* ( HoldRealNorm# “ 2  +  HoldlBNorm#"2 ) 


for  uX  s  0  to  fX  *  range  of  kasin( theta) 

u  =  uX/10 
incr  elementX 

SinSquareTheta#  =  ( u/ ( k#*a# ) ) " 2 


‘w"  ■N.’  ■-*  %’  ’  *  *  '  *.*  •• 
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'f  /  t  ■/  J 


Theta# 


atn( sqr( SinSquareTheta#/ ( 1 -SinSquareTheta# ) ) ) 


R#  =  z#/ ( f f min#*cos { The ta# ) )  '  calculate  relative  range 

if  elementX  =  0  then  g#  =  0  else  g#  =  z#*tan ( Theta# ) 

s#=( (l+cos( Theta# ) ) /R# ) "2_ 

*{RealOff#( elementX ) " 2  +  ImOf  f  # ( element% ) ' 2 ) / norm# 

if  elementX  =  0  then  save#  =  s# 


Iprint  using  "  #####.#  ";u; 


if  Real0ff#(0)  =  0  or  ImOff#(0> 
Iprint  "  **♦**»*»  ; 

elseif  RealOf f # ( e lementX ) =  0  or 
Iprint  "  ♦*****♦♦  ; 

else 

Iprint  using  "  #######.♦### 

end  i  f 


=  0  then 

ImOf f #( elementX  )  =  0  then 


;10*logl0(s#/8ave#); 


Iprint  using 


#.######  '’;g#: 


if  RealOf f #( elementX  )  =  0  or  ImOf f# ( elementX  )  =  0  then 
Iprint  " 
else 

Iprint  using  "  #.######  " ; s#* f f pwrden# 

end  i  f 


next  uX 

ParameterListOf f : 

1  print 

Iprint 

n 

H  parameter 

^  n  ^ 

Iprint 

using  " 

#.####" iH# 

Iprint 

«< 

sidelobe  ratio  ="; 

Iprint 

using  '' 

##.##♦  ":slr 

t 

Iprint 

Iprint 

"dB" 

M 

frequency  =” 

f 

Iprint 

using  " 

#.##### - 

" ; f  req#  ; 

Ipi.  int 
Iprint 

"GHz" 

#» 

diameter  = " ; 

Iprint 

using  " 

#.##### - 

" ; diam# ; 

Iprint 

"meters 

If 

Iprint 

If 

gain  =": 

Iprint 

using  " 

##  .  #### 

" ; gain; 

Iprint 

Iprint 

"dB" 

ft 

transmitted 

power  = 

Iprint 

using  " 

#.##♦### 

" ; power 

Iprint 

Iprint 

"kW" 

If 

off-axis  starting  po 

Iprint 

using  " 

#.#######'■'" 

~  ";z#; 

Iprint 

"meters 

If 

[♦A*] 


Iprint  "  error  Unit  ="; 

Iprint  using  "##.###"; 1 00*ErrLimOf f ; 
Iprint  "X" 
goto  Recycles 

DiscWriteOf  f : 


print 

print  "Job  complete; 
print  B$ 

open  B$  for  output  as 
print  #2, 
print  #2, 

print  #2,  :  print  #2, 

print  #2," 
print  #2," 
print  #2," 
print  #2,"kasin("; 
print  #2,  chr${233); 
print  #2 ,  " ) 
print  #2, 


writing  data  to  file  ' 

«2 


Off-Axis  Power  Density 


Relative 

Power  Distance 

Density  Off-Axis 


(decibels)  (meters) 


Data” 


Absolute" 

Power” 

Density" 


(mW/cm'‘2  )" 


elementX  =  -1  ’  normalize  wrt  ffmin# 

norm#  =  4* ( HoldRealNorm# ^ 2  +  Holdl mNorm# *  2 ) 


for  uX  =  0  to  fX  ’  range  of  kasin( theta) 

u  s  uX/10 
incr  elementX 

SinSquareTheta#  =  ( u/ ( k##a# ) ) ' 2 

Theta#  =  atn ( sqr ( S inSquareTheta# / ( 1 -S inSquareTheta# ) ) ) 

R#  =  z#/ ( f fmin#*cos ( Theta# ) )  ’  calculate  relative  range 

if  elementX  =  0  then  g#  =  0  else  g#  =  z#*tan( Theta# ) 

s#=( ( l+cos(Theta# ) )/R#)'2_ 

*(RealOff#(  elementX  )  2  +  ImOf  f  #  (  elementX  )  *  2  )  /  norm# 


if  elementX  =  0  then  save#  =  s# 

print  #2,  using  "  #####,#  ";u; 

if  RealOff#(0)  =  0  or  Im0ff#(0)  =  0  then 
print  #2,  "  ♦♦**♦*♦*  ; 

elseif  RealOff#(elementX)  =0  or  ImOf f# ( elementX )  =  0  then 
print  #2,  "  *♦*♦♦♦♦♦  ••  j 

else 

print  #2,  using  "#######.####  " ; 10*iogl0( s#/save#  )  ; 

end  if 


":g#; 


print  #2,  using 


#.f###«#^ 


if  RealOf f # ( elementX )  =  0  or  ImOf f # ( elementX )  =  0  then 
print  #2,  ’’ 
else 

print  #2,  using  "  ~ ~ s#* f f pwrden# 

end  i  f 

next  uX 

f . ParameterLi stOf f : 


print 

#2, 

print 

#2, 

tt 

H  parameter  ="  ; 

pr  i  nt 

#2, 

using 

tl 

« .««««" ;H« 

print 

♦  2, 

(1 

sidelobe  ratio  ="; 

print 

#2. 

using 

ft 

##.###  ";8lr; 

print 

#2, 

"dB" 

print 

#2, 

tt 

frequency  =" ; 

print 

#2, 

using 

fl 

#  .  ##### . ;  freq#  ; 

print 

#2, 

"GHz" 

print 

#2, 

diameter  ="  ; 

print 

#2, 

us  i  ng 

ft 

«  .  ##### . .  ;diam#  ; 

print 

♦  2, 

"meters 

ft 

print 

#2, 

M 

gain  =" : 

print 

#2, 

using 

If 

##.####  ";gain; 

print 

#2, 

"dB" 

print 

#2, 

'* 

transmitted  power  ="; 

print 

#2, 

using 

tt 

#.######  power; 

print 

#2, 

"kW" 

print 

#2, 

ft 

off-axis  starting  point 

print 

#2, 

using 

#.#######  '  ";z#; 

print 

♦  2, 

"meters 

print 

#2, 

" 

error  limit  ; 

print 

#2, 

using 

♦#.###": 100*ErrLimOf f ; 

print 

#2, 

"X" 

close 

#2 

Ref  ycleS  : 
print 

print  "Do  you  need  more  off-axis  data 

input  "using  the  current  parameters  (y/n)?  ",ans$ 

if  ans$  =  "y"  then 
goto  OffAxisStart 
elseif  ans$  =  "n"  then 
exit  if 
else 

goto  Recycles 
end  i  f 

Recycles : 


print 

input  "New  problem  (y/n)?  ",ans$ 


y"  then 


if  ans$  = 
els 
goto  Begin 

elseif  ans$  =  "n"  then 
exit  if 
else 

goto  Recycles 
end  i  f 
end 

def  fnF#{x#)  ’  used  in  Newton-Raphson  formula 

local  iX,w#,j#  ’  for  finding  the  H  parameter 

w#  =  0  :  j#  =  1E06 
for  iX  =  1  to  170 

w#=w#  +  X#*  (  2*iX)/(  2'' (  2*iX+l  )*fnFac#(  iX  )*fnFac#(  iX+1  )  ) 
if  abs(w#-j#)  <  tolerance  then  exit  for  else  j#  =  w# 
next  iX 

fnF#  =  w#  +  0.5  -  il# 
end  def 

def  fnG#(x#)  ’  returns  derivative  of  fnF#(x#) 

local  iX , y# , j# 
y#  =  0  ;  j#  =  1E06 
for  iX  =  1  to  170 

y#  =  y#  +  2*iX*x#' ( 2*iX-l )/( 2' ( 2*iX+l  )_ 

♦  fnFac#( iX)*fnFac#( iX+1  )  ) 

if  abs(y#-j#)  <  tolerance  then  exit  for  else  j#  =  y# 
next  iX 
fnG#  =  y# 
end  def 

def  fnFac#(xX)  ’  returns  the  factorial  of  a  number 

local  iXie# 
e*  =  1 

for  iX  =  1  to  xX 
e#  =  e#  *  iX 
next  iX 
fnFac#  =  e# 
end  def 

’  returns  sum  of  midpoint  evaluations 
def  f nMidSum# ( xX  )  ’  when  using  Trapezoidal  Rule  during 

local  inct , p# , iX , b# , d#  *  Romberg  integration 

inc#  =  1/xX  :  b#  =  0 
for  iX  =  1  to  (xX-1) 
p#  =  l-iX*inc# 
if  flagX  =  0  then 

d#  =  COB ( q#*p# " 2  ) 
else 

d#  =  8 i n ( q#  *p# " 2 ) 
end  i  f 

b#  =  b#  +  d#*p#*fnI#{pi#*H#*8qr(l-p#'2 ) )*fnJ#(u»p# ) 
next  iX 
fnMidSum#  =  b# 
end  def 
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def  fnl#(x#)  ’  returns  modified  Bessel  function 

local  iX,v#,j#  ’  of  the  first  kind,  order  zero 

V#  =  0  :  j#  =  1E06 
for  iX  =  1  to  170 

V#  =  V#  +  (x#/2)"(2*iX)/fnFac#( iX)'2 

if  abs(v#-j#)  <  tolerance  then  exit  for  else  j#  =  v# 
next  iX 
fnl#  =  V#  +  1 
end  def 

’  returns  Bessel  function  of  the  first  kind,  order  zero 
def  f  nJ# { X#  ) 

if  X#  >  25  then  fnJ#  =  fnA#(x#)_ 

♦  s in ( X# )  +  f nB# ( X#  )  ♦cos ( X#  )  :  exit  def 

local  iX,c#,m4 
c#  =  0  :  m#  =  1E06 
for  iX  =  1  to  170 

c#  =  c#  +  (-1  )''iX^(x#/2)*(2*iX)/fnFac#(  iX)"2 
if  abs(c#-m#)  <  tolerance  then  exit  for  else  m#  =  c# 
next  iX 
fnJ#  =  c#  +  1 
end  def 

’  fnA# , fnB# , fnP# , fnQ#  used  when  x#  >  25  for  fnJ#(x#) 
def  fnA#(x#)  =  (fnP#(x#)  -  f nQ# ( x#  )  ) /sqr ( pi#*x#  ) 
def  fnB#(x#)  =  (fnP#(x#)  +  fnQ# ( x# ) ) /sqr ( pi#*x# ) 
def  f nP# ( X#  ) 

local  wX , cX , check# , mX , sum# , iX , product# 

mX  =  0  wX  =  -1  :  cX  =  0  :  check#  =  1E06  :  sura#  =  0 

Rout ineP : 

w%  =  wX  +  4  :  cX  =  cX  +  2  :  mX  =  mX  +  1  :  product#  =  1 
for  iX  =  1  to  wX  step  2 

product#  =  product#  *  iX"2 
next  iX 

product#  =  product#*  (- 1  ) ''mX 
product#  =  produc t# / f nFac# { cX ) 
product#  =  product#/ {  8*x#  ) ''cX 
sura#  =  sum#  +  product# 

if  abs(sum#  -  check#)  <  tolerance  then 
sum#  =  sum#  -f  1 
else 

check#  =  sum# 
goto  RoutineP 
end  i  f 
fnP#  =  sum# 
end  def 
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def  fnQ#(x#) 

local  wX , cX , check# , mX , sum# , iX , product# 

mX  =  0  :  wX  =  -3  cX  =  -1  :  check#  =  1E06  sum#  =  0 

Rout ineQ : 

wX  =  wX  +  4  cX  =  cX  +  2  mX  =  mX  +  1  product#  =  1 
for  iX  =  1  to  wX  step  2 

product#  =  product#  •  iX~2 
next  iX 

product#  =  product#* { -1 ) “mX 
product#  =  product#/fnFac# ( cX ) 
product#  =  product#/ ( 8*x# ) "c% 
sum#  =  sum#  +  product# 

if  abs(sum#  -  check#)  <  tolerance  then 
exit  if 
else 

check#  =  sum# 
goto  RoutineQ 
end  i  f 
fnQ#  =  sum# 
end  def 


********** 


GLOSSARY 

sir  .  antenna  sidelobe  ratio 


H  parameter  for  one-parameter 
distribution  developed  by  Hansen 


pi#  .  3.1416 

est#  . .  iterative  value  of  pi#*H*  used  in  loop 

employing  Newton’s  method  to  solve  for 
the  H  parameter 

tolerance  .  accuracy  to  which  internal  calculations 

are  computed 

freq#  .  operating  frequency 

lambda#  .  wavelength 


propagation  constant 


diam#  .  antenna  diameter 

RangeMin#  .  minimum  axial  range 

ffmin#  .  minimum  far  field  distance 


a#  .  antenna  radius 

gain  .  antenna  power  gain  in  decibels 

power  .  transmitted  power  (peak  or  average 

depending  on  the  hazard  criteria) 

ffpwrden#  .  power  density  at  ffmin# 

ErrLimOn  .  maximum  error  accepted  by  Romberg 

integration  routine  when  calculating 
axial  power  density 

RealOn#  .  linear  array  element;  stores  real  part 

of  integral  during  axial  computations 

ImOn#  .  linear  array  element;  stores  imaginary 

part  of  integral  during  axial 
computations 

u  .  k#*a#*sin( theta ) 

R,R#  .  relative  distance  from  aperture  origin 

to  field  point;  normalized  with  respect 
to  ffmin# 

Romberg#  .  12x12  array  element;  stores  all  elements 

of  the  Romberg  matrix  until  the  error 
limit  is  reached 

z#  .  axial  range  at  which  off-axis 

computations  begin 

ThetaMax#  .  maximum  off-axis  angle 

Theta#  .  present  off-axis  angle 

f  .  value  of  k#*a#*sin( Theta#  )  to  proceed  to 

off-axis 

ErrLimOff  .  maximum  error  accepted  by  Romberg 

integration  routine  when  calculating 
off-axis  power  density 

RealOff#  .  linear  array  element;  store  real  part  of 

integral  during  off-axis  computations 

ImOff#  .  linear  array  element;  store  imaginary 

part  of  integral  during  off-axis 
computations 
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Appendix  B:  POWERDEN  Axial  Output  Data 


Axial  Power  Density  Data 


Relative 

Absolute 

Power 

Absolute 

Power 

Relative 

Density 

Range 

Density 

Range 

( decibels ) 

( aeters ) 

( mW/cm' 2 ) 

1.0000 

0.0000 

2 .400000E+03 

1 .OOOOOOE+01 

0.9900 

0.0862 

2.376000E+03 

1 .020037E+01 

0.9800 

0.1732 

2.352000E+03 

1 .040680E+01 

0.9700 

0.2611 

2.328000E+03 

1 .061952E+01 

0.9600 

0.3498 

2.304000E+03 

1 .083881E+01 

0.9500 

0.4395 

2.280000E+03 

1 . 106492E+01 

0.9400 

0.5301 

2.256000E+03 

1 . 129814E+01 

0.9300 

0.6216 

2.232000E+03 

1 . 153878E+01 

0.9200 

0.7141 

2.208000E+03 

1 . 178715E+01 

0.9100 

0.8076 

2. 184000E+03 

1 .204357E+01 

0.9000 

0.9020 

2.160000E+03 

1 .230842E+01 

0.8900 

0.9975 

2. 136000E+03 

1 .258204E+01 

0.8800 

1.0940 

2. 112000E+03 

1.286483E+01 

0.8700 

1.1916 

2.088000E+03 

1 . 315721E+01 

0.8600 

1.2903 

2.064000Et03 

1 .345960E+01 

0.8500 

1.3901 

2.040000E+03 

1.377248E+01 

0.8400 

1.4911 

2.016000E+03 

1 .409632E+01 

0.8300 

1.5932 

1.992000E+03 

1 .443164E+01 

0.8200 

1.6964 

1.968000E+03 

1 .477898E+01 

0.8100 

1.8010 

1 .944000E+03 

1 . 513893E+01 

0.8000 

1.9067 

1 . 920000E+03 

1 . 551209E+01 

0.7900 

2.0137 

1 .896000E+03 

1 .589912E+01 

0.7800 

2.1221 

1 . 872000E+03 

1 .630071E+01 

0.7700 

2.2317 

1 .848000E+03 

1 .671758E+01 

0.7600 

2.3428 

1 .824000E+03 

1 . 715052E+01 

0.7500 

2.4552 

1 .800000E+03 

1 . 760035E+01 

0.7400 

2.5691 

1 . 776000E+03 

1 .806796E+01 

0.7300 

2.6844 

1.752000E+03 

1 .855429E+01 

0.7200 

2.8013 

1 . 728000E+03 

1 .906034E+01 

0.7100 

2.9197 

1.704000E+03 

1.958719E+01 

0.7000 

3.0397 

1.680000E+03 

2.013597E+01 

0.6900 

3.1614 

1 .656000E+03 

2.070790E+01 

0.6800 

3.2847 

1 .632000E+03 

2.130430E+01 

0.6700 

3.4097 

1 .608000E+03 

2 . 192656E+01 

0.6600 

3.5365 

1 .584000E+03 

2.257618E+01 

0.6500 

3.6651 

1 . 560000E+03 

2.325477E+01 

0.6400 

3.7956 

1 .636000E+03 

2 . 396406E-»-01 

0.6300 

3.9280 

1.512000E+03 

2.470589E-t-01 

0.6200 

4.0624 

1.488000E-t-03 

2.548226E+01 

0.6100 

4.1988 

1.464000E+03 

2.629532E-»^01 

0.6000 

4.3373 

1 .440000E+03 

2. 714739E+01 

0.5900 

4.4779 

1.416000E+03 

2.804096E-f01 
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0.5800 

4.6208 

1.392000E+03 

2.897871E+01 

0.5700 

4 . 7659 

1 . 368000E+03 

2.996357E+01 

0.5600 

4.9134 

1.344000E+03 

3 .099869E+01 

0.5500 

5.0634 

1 . 320000E+03 

3.208747E+01 

0.5400 

5.2158 

1 .296000E+03 

3.323363E+01 

0.5300 

5.3708 

1 . 272000E+03 

3.444118E+01 

0.5200 

5.5284 

1.248000E+03 

3. 571449E+01 

0.5100 

5.6889 

1 . 224000E+03 

3. 705830E+01 

0.5000 

5.8521 

1 .200000E+03 

3.847780E+01 

0.4900 

6.0183 

1 . 176000E+03 

3.997863E+01 

0.4800 

6.1875 

1 . 152000E+03 

4. 156697E+01 

0.4700 

6.3598 

1  . 128000E+03 

4. 324954E+01 

0.4600 

6.5354 

1 . 104000E+03 

4.503373E+01 

0.4500 

6.7143 

1 .080000E+03 

4.692762E+01 

0.4400 

6.8966 

1 .056000E+03 

4.894008E+01 

0.4300 

7.0826 

1 .032000E+03 

5. 108086E+01 

0.4200 

7.2722 

1 . 008000E+03 

5.336068E+01 

0.4100 

7.4657 

9.840000E+02 

5.579131E+01 

0.4000 

7.6631 

9.600000E+02 

5 .838575E+01 

0.3900 

7.8646 

9.360000E+02 

6. 115836E+01 

0 . 3800 

8 .0703 

9 . 120000E+02 

6.412492E+01 

0.3700 

8.2803 

8. 880000E+02 

6. 730292E+01 

0.3600 

8.4949 

8.640000E-t-02 

7.071166E+01 

0.3500 

8.7141 

8.400000E+02 

7.437247E+01 

0.3400 

8.9381 

8. 160000E+02 

7.830890E+01 

0.3300 

9.1670 

7.920000E+02 

8.254696E+01 

0.3200 

9.4009 

7.680000E+02 

8. 711535E+01 

0.3100 

9.6400 

7.440000E+02 

9.204558E+01 

0.3000 

9.8844 

7.200000E+02 

9. 737226E+01 

0.2900 

10.1340 

6.960000E+02 

1 .031331E+02 

0.2800 

10.3889 

6. 720000E+02 

1 .093690E+02 

0.2700 

10.6492 

6.480000E+02 

1 . 161238E+02 

0.2600 

10.9147 

6.240000E+02 

1 . 234437E+02 

0.2500 

11  .  1852 

6.000000E+02 

1 . 313763E+02 

0.2400 

11.4601 

5. 760000E+02 

1.399627E+02 

0.2300 

11 . 7394 

5.520000E+02 

1 .492590E+02 

0.2200 

12.0220 

5.280000E+02 

1 .592955E+02 

0.2100 

12.3069 

5.040000E+02 

1 . 700928E+02 

0.2000 

12.5922 

4.800000E+02 

1 .816440E+02 

0.1900 

12.8757 

4.560000E+02 

1 .938961E+02 

0.1800 

13.1539 

4.320000E+02 

2.067216E+02 

0.1700 

13.4217 

4.080000E+02 

2. 198743E+02 

0.1600 

13.6721 

3.840000E+02 

2.329211E+02 

0.1500 

13.8941 

3.600000E+02 

2.451403E+02 

0.1400 

14.0717 

3.360000E+02 

2.553681E+02 

0.1300 

14.1794 

3. 120000E+02 

2.617814E+02 

0.1200 

14.1766 

2.880000E+02 

2.616136E+02 

0.1100 

13.9944 

2.640000E+02 

2.508631E+02 

0.1000 

13.5074 

2.400000E+02 

2.242541E+02 

0.0990 

13.4340 

2.376000E+02 

2.204970E+02 

0.0980 

13.3550 

2.352000E+02 

2. 165187E+02 

0.0970 

13.2698 

2.328000E+02 

2. 123147E+02 

0.0960 

13.1781 

2.304000E+02 

2.078806E+02 

0.0950 

13.0795 

2. 280000E+02 

2.032129E-»^02 

b 

-I'"' 

/v*, 

m 

i  * 

k  0.0940 

12.9734 

2.256000E+02 

1 .983086E+02 

•  *  1 

m.  ■ 

b  0.0930 

12.8595 

2. 232000E+02 

1 .931761E+02 

f  0.0920 

12.7368 

2.208000E+02 

1 .877923E+02 

L  0.0910 

12.6047 

2. 184000E+02 

1 .821674E+02 

•■4 

1  0.0900 

12.4626 

2. 160000E+02 

1.763022E+02 

;•  0.0890 

12.3096 

2. 136000E+02 

1 . 701984E+02 

y  0.0880 

12.1447 

2,112000E+02 

1 .638596E+02 

i  0.0870 

11.9670 

2.088000E+02 

1 . 572906E+02 

K  0.0860 

11 . 7753 

2,064000E+02 

1 .504985E+02 

<•,  . 

2  0.0850 

11 . 5683 

2.040000E+02 

1 .434921E+02 

fP  0.0840 

11.3444 

2. 016000E+02 

1 . 362829E+02 

r-  % 

0.0830 

11.1020 

1 .992000E+02 

1 ,288847E+02 

‘V 

r’  0.0820 

10.8391 

1 . 968000E+02 

1 .213145E+02 

*  • 

V  * 

^  0.0810 

10.5535 

1 . 944000E+02 

1 . 135922E+02 

f-  0.0800 

10.2425 

1 . 920000E+02 

1 .057416E+02 

i  0.0790 

9.9030 

1 . 896000E+02 

9. 779032E+01 

1  0.0780 

9.5313 

1 .872000E+02 

8. 976996E+01 

V 

0.0  7  70 

9.1231 

1 .848000E+02 

8. 171710E+01 

y  0.0760 

8.6731 

1 .824000E+02 

7 . 367328E+01 

0.0750 

8.1747 

1 . 800000E+02 

6.568555E+01 

>•  0.0740 

7.6198 

1 . 776000E+02 

5 . 780662E+01 

‘•X 

0.0730 

6.9980 

1 . 752000E+02 

5.009560E+01 

s 

•  0.0720 

6.2959 

1 . 728000E+02 

4.261813E+01 

y  0.0710 

5.4958 

1 . 704000E+02 

3.544680E+01 

y 

^  0.0700 

4.5729 

1 .680000E+02 

2.866114E+01 

0.0690 

3.4924 

1 .656000E+02 

2.234827E+01 

••x 

C:  0.0680 

2.2017 

1 .632000E  +  02 

1 .660233E+01 

’-■V 

i  0.0670 

0.6165 

1 .608000E+02 

1 . 152522E+01 

.  0.0660 

-1.4115 

1 .584000E+02 

7.225278E+00 

*  41  ’ 

0.0650 

-4 .1823 

1.560000E+02 

3.817386E+00 

0.0640 

-8.4713 

1 .536000E+02 

1 .421895E+00 

'  0.0630 

-17.8700 

1.512000E+02 

1 .633047E-01 

s*  » 

^  0.0620 

-17 . 7311 

1 .488000E+02 

1 .686141E-01 

j  0.0610 

-8.0551 

1 .464000E+02 

1.564900E+00 

P  0.0600 

-3.4908 

1 .440000E+02 

4.476353E+00 

. 

0.0590 

-0.4477 

1 .416000E+02 

9.020411E+00 

^  0.0580 

1,8476 

1 .392000E+02 

1 .530248E+01 

0.0570 

3.6944 

1 . 368000E+02 

2.341234E+01 

:•  0.0560 

5.2393 

1 .344000E+02 

3.341413E+01 

^  0.0550 

6.5648 

1.320000E+02 

4 . 534000E+01 

■  0.0540 

7.7218 

1 . 296000E+02 

5.918129E+01 

►:  0.0530 

8.7432 

1 . 272000E+02 

7.487274E+01 

P  0.0520 

9.6514 

1.248000E+02 

9.228749E+01 

f  0.0510 

10.4618 

1 .224000E+02 

1 . 112195E+02 

•:  0.0500 

11 . 1850 

1 .200000E+02 

1 .313707E+02 

^  0.0490 

11.8281 

1 . 176000E+02 

1 . 523375E+02 

s*  s 

J  0.0480 

12.3955 

1 . 152000E+02 

1  .  735984E+02 

.s 

V  0.0470 

12.8893 

1 . 128000E+02 

1 . 945051E+02 

1;  0.0460 

13.3098 

1 . 104000E+02 

2 . 142772E+02 

'•  0.0450 

13.6547 

1 .080000E+02 

2.319897E+02 

0.0440 

13.9209 

1 .056000E+02 

2.466537E+02 

/  'j 
/•j 

•  0.0430 

14.1022 

1.032000E+02 

2. 571714E+02 

•  0.0420 

14 . 1892 

1 .008000E+02 

2.623741E+02 

o 

o 

o 

14.1698 
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9.840000E+01 

2 .612019E-t-02 

s' 

\  .s 
S'.N 

^  i. «-  fc.v ■  MVaViiVyHvyi4y_>  r-» 

'i'  S' 

1 

0.0400 

14.0268 

9.600000E+01 

2.527415E+02 

0.0390 

13.7364 

9. 360000E+01 

2.363942E+02 

0.0380 

13.2646 

9. 120000E+01 

2. 120604E+02 

0.0370 

12.5612 

8.880000E+01 

1 .803533E+02 

0.0360 

11.5466 

8. 640000E+01 

1 .427769E+02 

0.0350 

10.0919 

8.400000E+01 

1 .021381E+02 

0.0340 

7.9388 

8. 160000E+01 

6.221326E+01 

0.0330 

4.4881 

7. 920000E+01 

2.810652E+01 

0.0320 

-2.4743 

7.680000E+01 

5.656799E+00 

0.0310 

-11 . 7133 

7.440000E+01 

6. 740233E-01 

0.0300 

2.4552 

7.200000E+01 

1 . 760035E+01 

0.0290 

7.6077 

6. 960000E+01 

5 . 764668E+01 

0.0280 

10.6693 

6. 720000E+01 

1 . 166618E+02 

0.0270 

12.6341 

6.480000E+01 

1 .834039E+02 

0.0260 

13.7928 

6.240000E+01 

2 . 394866E+02 

0.0250 

14.1953 

6.000000E+01 

2.627413E+02 

0.0240 

13.7216 

5. 760000E+01 

2.355940E+02 

0.0230 

11 . 9893 

5.520000E+01 

1 . 580988E-I-02 

0.0220 

7.8065 

5.280000E+01 

6.034628E+01 

0.0210 

-8.3348 

5.040000E+01 

1 .467297E+00 

0.0200 

5.8521 

4.800000E+01 

3.847767E+01 

0.0196 

9.1419 

4.704000E+01 

8.207168E+01 

0.0192 

11 . 3248 

4.608000E+01 

1 .356690E+02 

0.0188 

12.8017 

4 .512000E+01 

1 . 906204E+02 

0.0184 

13.7290 

4.416000E+01 

2.359914E+02 

0.0180 

14.1622 

4 . 320000E+01 

2 . 607458E+02 

0.0176 

14.0827 

4.224000E+01 

2.560185E+02 

0.0172 

13.4052 

4. 128000E+01 

2.190368E+02 

0.0168 

11.9255 

4.032000E+01 

1 .557942E+02 

0.0164 

9.1520 

3.936000E+01 

8.226215E+01 

0.0160 

3.4518 

3.840000E+01 

2.214004E+01 

0.0156 

-19.7247 

3 . 744000E+01 

1 .065454E-01 

0.0152 

4.9309 

3. 648000E+01 

3.112377E+01 

0.0148 

10.3810 

3.552000E+01 

1 .091689E+02 

0.0144 

13.0526 

3.456000E+01 

2.019570E+02 

0.0140 

14 . 1405 

3.360000E+01 

2 . 594476E+02 

0.0136 

13.7976 

3.264000E+01 

2.397512E+02 

0.0132 

11.5793 

3. 168000E+01 

1 .438582E+02 

0.0128 

5.3188 

3.072000E+01 

3.403135E+01 

0.0124 

-3.7740 

2.976000E+01 

4. 193701E+00 

0.0120 

9.8842 

2. 880000E+01 

9. 736940E+01 

0.0116 

13.6454 

2. 784000E+01 

2.314964E+02 

0.0112 

13.9155 

2.688000E+01 

2.463492E-«-02 

0.0108 

10.0464 

2.592000E+01 

1 .010746E+02 

0.0104 

-16.2037 

2.496000E+01 

2.396779E-01 

0.0100 

11 . 1850 

2 . 400000E+01 

1 .313708E+02 

H  parameter 

=  0.0000 

sidelobe  ratio  =  17.570 

dB 

frequency  = 

1 .OOOOOE+01 

GHz 

diameter  =  6 

.OOOOOE+00 

meters 

gain  =  55.5860E+00  dB 

transmitted 

power  =  2.000000E+01  kW 

error  limit 

=  O.IOOX 
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.  vJ-y-V-V- 

Off-Axis  Power  Density  Data 


kasin( 6  ) 

Relative 
Power 
Density 
( decibels ) 

Distance 

Of  f-Axis 
( meters ) 

0.0 

0.0000 

O.OOOOOOE+00 

0.1 

-0.0109 

1 .909859E-01 

0.2 

-0.0435 

3.819719E-01 

0.3 

-0.0979 

5. 729579E-01 

0.4 

-0.1743 

7.639439E-01 

0.5 

-0.2728 

9.549300E-01 

0.6 

-0.3937 

1  . 145916E+00 

0.7 

-0.5373 

1 . 336902E  +  00 

0.8 

-0.7040 

1.527889E+00 

0.9 

-0.8941 

1 . 718875E+00 

1.0 

-1 . 1083 

1 .909862E+00 

1  .  1 

-1 . 3470 

2. 100849E+00 

1.2 

-1.6110 

2.291835E+00 

1.3 

-1.9010 

2.482822E+00 

1 .4 

-2.2180 

2.673810E+00 

1.5 

-2.5629 

2.864797E+00 

1.6 

-2.9368 

3.055785E+00 

1 . 7 

-3.3412 

3.246773E+00 

1.8 

-3.7773 

3.437761E+00 

1 . 9 

-4.2470 

3.628749E+00 

2.0 

-4.7522 

3.819738E+00 

2 . 1 

-5.2950 

4.010727E+00 

2.2 

-5.8779 

4.201716E+00 

2.3 

-6.5040 

4 . 392706E+00 

2.4 

-7.1764 

4.583696E+00 

2.5 

-7.8989 

4.774686E+00 

2.6 

-8.6759 

4.965677E+00 

2.7 

-9.5124 

5. 156668E+00 

2.8 

-10.4136 

5 . 347659E+00 

2.9 

-11 . 3855 

5.538651E+00 

3.0 

-12.4336 

5.729643E+00 

3.1 

-13.5625 

5.920636E+00 

3.2 

-14 .7730 

6. 111629E+00 

3.3 

-16.0571 

6.302623E+00 

3.4 

-17 . 3880 

6.493617E+00 

3.5 

-18.7035 

6.684611E+00 

3.6 

-19.8855 

6.875606E+00 

3.7 

-20.7619 

7.066602E+00 

3.8 

-21 .1745 

7.257598E+00 

3.9 

-21.0953 

7.448595E+00 

4.0 

-20.6512 

7.639592E+00 

4.1 

-20.0227 

7.830590E+00 

Absolute 
Power 
Density 
( mW/cm" 2  ) 

3.847780E+01 
3.838170E+01 
3.809463E+01 
3. 762018E+01 
3.696430E+01 
3.613515E+01 
3 . 514299E+01 
3 . 399998E+01 
3 . 271997E+01 
3. 131830E+01 
2.981148E+01 
2 . 821696E+01 
2 . 655280E+01 
2.483743E+01 
2 . 308931E+01 
2 . 132665E+01 
1 .956714E+01 
1 . 782768E+01 
1 .612422E+01 
1 .447135E+01 
1 . 288232E+01 
1 . 136878E+01 
9. 940688E+00 
8.606241E+00 
7.371804E+00 
6.241969E+00 
5. 219358E+00 
4 . 304985E+00 
3.498226E+00 
2.796808E+00 
2. 197087E+00 
1 .694163Et00 
1 .282057E+00 
9.538995E-01 
7 .021168E-01 
5. 186322E-01 
3.950542E-01 
3.228655E-01 
2.936003E-01 
2.990069E-01 
3.311966E-01 
3.827723E-01 
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4 . 2 

-19.3496 

8.021588E+00 

4.469383E-01 

» 

* « * 
•'.V 

4.3 

-18.7123 

8.212588E+00 

5. 175875E-01 

4.4 

-18.1483 

8 .403587E+00 

5.893656E-01 

4.5 

-17.6717 

8. 594587E+00 

6.577157E-01 

4.6 

-17.2854 

8.785588E+00 

7. 188985E-01 

1 

* » • 

4 . 7 

-16.9872 

8.976590E+00 

7.699923E-01 

4.8 

-16.7733 

9. 167593E+00 

8.088738E-01 

4.9 

-16.6395 

9. 358595E+00 

8. 341803E-01 

v: 

5.0 

-16.5813 

9.549599E+00 

8.454284E-01 

5.1 

-16.5985 

9. 740603E+00 

8.420955E-01 

>A 

5.2 

-16.6863 

9.931608E+00 

8.252479E-01 

\ 

5.3 

-16.8443 

1 .012261E+01 

7.957519E-01 

5.4 

-17.0724 

1 .031362E+01 

7,550382E-01 

5.5 

-17.3712 

1 .050463E+01 

7.048384E-01 

*  *  < 

5.6 

-17.7424 

1 .069564E+01 

6.470950E-01 

5.7 

-18.1889 

1 .088665E+01 

5.838726E-01 

. ' 

5.8 

-18.7149 

1 . 107766E+01 

5. 172739E-01 

5.9 

-19.3261 

1 . 126867E+01 

4 .493645E-01 

6.0 

-20.0303 

1 . 145968E+01 

3.821040E-01 

6.1 

-20.8376 

1 . 165069E+01 

3 . 172867E-01 

>■; 

6.2 

-21 . 7614 

1  .  184170E+01 

2.564899E-01 

» » «• 

6.3 

-22.8189 

1 .203272E  +  01 

2.010560E-01 

6.4 

-24.0326 

1 .222373E+01 

1 . 520388E-01 

1  \ 

6 . 5 

-25.4298 

1 .241475E+01 

1 . 102111E-01 

•* 

.‘/w- 

6.6 

-27.0416 

1 .260577E+01 

7.604226E-02 

. 

6.7 

-28.8861 

1 .279678E+01 

4.972785E-02 

6.8 

-30.9121 

1 .298780E+01 

3. 118907E-02 

6.9 

-32.8203 

1 .317882E+01 

2.009934E-02 

•V-4 

7.0 

-33.8331 

1.336984E+01 

1 . 591848E-02 

1 

•  VJ 

7.1 

-33.3188 

1 .356087E+01 
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sidelobe  ratio  =  17.570  dB 

frequency  =  l.OOOOOE+01  GHz 

diameter  =  6.00000E+00  meters 

gain  =  55.5860E+00  dB 

transmitted  power  =  2.000000E+01  kW 

off-axis  starting  point  =  1 . 2000000E+03  meters 

error  limit  =  O.IOOX 
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Appendix  D:  Derivation  of 


The  scalar  Kirchoff  diffraction  integral  with  the  FSA 
approximations  is 


f  ra  -jkp2 

Ep  =  B(6)1  /G(p)exp - exp[  jkpsin  (  ^  )  cos  ( 0 -0  )  ]  pdpd<<>  (6) 

o  r» 


where 


B{^>)  =  (l+cos(tf))  (jk/47rR)  exp{-jkR) 


and  all  other  terms  are  as  defined  originally  in  Chapter  2. 
Since  only  relative  amplitude  is  desired,  the  constant  jk 
and  the  phase  term  can  be  dropped  from  B(S)  resulting  in 


lEpl~C(e)U  jG(p)exp 


exp[ jkpsini $ ) cos (0-0) ]pdpd0  ( 18  ) 


where 


c(  e )  =  [  i+cos(  e )  ]/47rR 


With  0  arbitrarily  set  to  zero  since  the  radiation  field  is 
circularly  symmetric,  the  integration  with  respect  to  0  can 
be  performed  through  use  of  the  identity  (Stutzman  and 
Thiele,  1981:416,566) 


Jo(x)  =  (1/2T)  /  expl jxcos ( a ) ]  da 


‘  _  n  -  •  «>  «  •  mV*.  »  *  *  ^  ' 


k  s 


so 
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at  i  ve 

field 

strength 

is 

1  1 

-  jkp2 

|ep|  ~  d(« 

) 

J  G(P 

)  exp 

— 

Jo 

[ kpsini 0 ) J  p  dp 

2R 

where 

D(  0  )  = 

[  1 

+  COS  (  0 

)  ]/2R 

If 

the  subst 
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X  =  P/a,  A  = 

R/ 

{2D2/X),  and 

u 

=  kas in (  0 ) 

ar 

e  made 

,  and 

the  CO 

nstants  appearing  ( 

(2 


the  integral  are  dropped  then 
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